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I. INTRODUCTION 


A. FRACTAL GEOMETRY AS A DESCRIPTION OF NATURE 

Fractal Geometry originated in the late 1970’s with Benoit 
Mandelbrot (ref. 1]. Mandelbrot felt that lines, circles and 
spheres of Euclidean geometry were unable to describe such 
patterns of nature as clouds, coastlines, trees and mountains. 
He developed the concept of fractal geometry to describe these 
patterns of nature. Although the concept of Fractal Geometry 
was a new concept it includes earlier work of Peano, Cantor 
and Koch [ref. 1].  Peano’s space filling curves, Cantor’s 
dust concept and the Koch curves have all been incorporated 
into the theory of fractal geometry. A Peano curve is a 
mapping of a two dimensional curve on the plane (2D) which 
completely fills the two dimensional space. The Cantor dust 
concept, also called Cantor’s middle thirds set is a 
topological set of dimension 0. The set is created on the 
closed interval by dividing the interval [0,1] into three 
pieces and then removing the middle third. Each of the 
remaining pieces are further divided and then the middle 
thirds removed in the same way and so on to infinity. Koch 
curves are particularly relevant to this study because the 
Koch snowflake is the first application of the midpoint 
displacement technique used here to generate fractal mountains 


Meetic, pp. 1-4.) 


B. KOCH SNOWFLAKE 

The Koch snowflake is one of the most well known fractal 
curves. The generation of the snowflake illustrates part of 
the approach taken in this study. The snowflake is created 
by starting with a basic Euclidean shape. A line segment, 
such as AB in Figure 1.la, is used as an initiator. The 
midpoint of the initiator (or base line segment), C’, is 
displaced to a point C as in Figure 1.1b. The initiator AB 
is replaced by the pair of lines AC and CB. The midpoints of 
line segments AC and CB are then displaced back to the 
original line at points D and F respectively. By connecting 
D to C by the line DC and C to F by the line CF the resulting 
polygonal line, ADCFB shown in Figure 1.1c, is called the 
generator. The process used to create the generator is called 


the midpoint displacement technique. 


To create the fractal or Koch snowflake, the process of 
creating generators out of each initiator is performed 
recursively. A process or a function is said to be recursive 
if, in the course of its execution, the function issues a call 
to itself [ref. 2, pp. 37-38]. In the case of the Koch 
snowflake, as in any recursive process, creating generators 
from each new initiator is, in theory, an infinite process. 
However, repeating a process or function infinitely often in 
this way is not possible, except theoretically. On finite 


equipment such as a computer there must be a criterion 





Figure 1.1 Koch Initiator and Generator 


developed to stop the process. The criterion for stopping is 
developed with the limitations of the computer being used in 
the recursion kept in mind. For example, some smallest line 


segment length may be incorporated into a stopping criterion. 


C. THREE DIMENSIONAL EXTENSION 

In this study, the midpoint technique has been extended 
to apply to three dimensional objects. To extend the midpoint 
displacement technique from two to three dimensions, three 
line segments in the form of a triangle are used as 
initiators. This triangle, on which the fractal mountain is 
created, is called a base triangle and, as can be seen in 
Figure 1.2, is placed in the xz plane. The right handed 


coordinate system shown in Figure 1.2 is the coordinate system 


that is used in this study. The y coordinate is viewed as the 
vertical direction and negative values of z are taken as going 
into the plane as the viewer perceives the vectors in the 
coordinate system. The base triangle can be taken to be any 
shape or size. By the use of many base triangles, the 
development of a fractal mountain scene rather than just a 
single mountain as in the case of one base triangle is 


possible. 





Figure 1.2 Base Triangle and Coordinate 


System 


De MOTIVATION 
Most of the earlier literature on fractal mountains has 


been developed from an artistic approach. The artistic 


approach gives the desired mountain views, but it is difficult 
for someone other than the artist to repeat the results as no 
systematic approach is available. More recent literature has 
included some algorithms on the different techniques for 
creating fractal mountains. In order to formalize the 
tinkering or artistic tricks used in the creation of a 
realistic picture, fractal mountains must be developed using 
interactive algorithms that can be easily understood and 
duplicated by others. Further improvements may then be added 


to enhance and expand the possible types of fractal mountains. 


E. SUMMARY OF THE REST OF THE STUDY 

In Chapter II, a literature review of the historical 
development of fractal mountains is given. The basic theory 
and algorithms used to create the fractals are discussed and 
a stopping criterion based on machine limitations on line 
segment length also appears there. In Chapter III, the main 
technical limitations to the quality of fractal mountain 
creation, granularity and gaps, are discussed. In Chapter IV, 
the problems connected with algorithms used in solving the 
granularity and gapping problems are discussed. In Chapter 
V, the final implementation used in this study to achieve a 
natural look in the fractal mountain generation is given. 
The final chapters explain the Interactive System for Fractal 
Mountains (ISFM). ISFM is implemented by pop-up menus which 


provide the options needed to change the texture of the 


fractal mountain and to change the light source position for 
either a mountain or a mountain range. The contour for the 
mountain range can also be defined interactively. The 
description of ISFM is followed by actual pictures showing 


fractal mountains. 


II. FRACTAL MOUNTAINS 


A. HISTORIC BACKGROUND 

The application of fractals has become more extensive 
since the fractal concept was originated in the late 1970’s. 
Mort La Brecque [ref. 3] cites applications ranging from 
chemical reactions, protein behavior, mechanical and 
electrical system behavior to natural objects. The natural 
objects created by fractals include galactic clusters, 
earthquake patterns, rainfall, clouds and landscapes. The 
results obtained in each of these applications have been 
enhanced by advances achieved in computer graphics. With the 
continued expanding technology of the computer, images created 
by fractals help to illustrate and explain many different 
aspects of all the applications discussed. 

Fractals were first developed using the ideas of self- 
Similarity [ref. 1,4,5,6] and self-affinity (ref. 5,6]. 
Michael F. Goodchild [ref. 4] defines self-similarity as a 
property of certain curves whereby each part of the curve is 
indistinguishable from the whole curve and by which, on each 
different level of scale, the finer structure of the curve 


resembles the grosser structure. Self-affinity is defined by 


Richard F. Voss. [ref. 5] Consider a set S = {x,y,.-..} of 
points defined, as x=(x;,....,X%¢), in a Euclidean space of 
dimension E. Under an affine transformation, each of the E 


coordinates of x may be scaled by a different ratio r = 
(Y1,Y2, +.» »Yeg). The set S is transformed to rS with points at 
r(x) = (Y4X1,++.+.,YeXe). A bounded set S is called self-affine 
when S is the union of N distinct (non-overlapping) subsets 
each of which is congruent to rS. Self-similarity and self- 
affinity are essentially geometric scaling properties. The 
Koch snowflake discussed in Chapter I is an example of a curve 
that has self-similarity. That is, as we iteratively produce 
the snowflake by repeated operations of the generator on each 
new initiator for each level of scale, the same structure for 
each level is repeated yet finer. 

There are several other effective fractal methods which 
can be used to model natural objects. Demko, Hodges and 
Naylor (ref. 7] produce iterated function systems for 
generating fractals in two dimensions (2D). These iterated 
systems can produce many of the classic 2D deterministic 
fractals found in Mandelbrot [ref. 1]. Another fractal method 
used to represent nature is based on the recursive subdivision 
method [ref. 8,9,10,11,12]. Amburn, Grant and Whitted [ref. 
8] indicate that the recursive subdivision procedure has three 
components, namely: 

1) framework 
2) subdivision equation 
3) termination criterion 

The framework is the overall organization of the procedure 

and controls the order of surface subdivisions. The 


subdivision equations describe the shape of the original 


polygon and the shapes of the iterated polygons at each level. 
The termination criterion is used in ending the fractal 
subdivision process. Some of the subdivision equations 
described by Gavin Miller [ref. 9] can result in triangle-edge 
subdivisions, diamond-square subdivisions or square-square 
subdivisions. Loren C. Carpenter [ref. 10] and Alan Ray Smith 
(ref. 11] describe yet another subdivision equation based on 
a triangular subdivision. The termination criterion is 
affected by the computer that is being used and is employed 
to stop the recursion process at the level in which the 
criterion is met. 

Carpenter (ref. 12] describes an additional subdivision 
method which divides a fractal region of finite size in such 
a way as to maintain self-similarity in the object created and 
thereby the statistical properties of the object. Carpenter’s 
method determines the midpoint of a fractal curve that is to 
be subdivided by a constrained random process. Fournier, 
Fussell and Carpenter [ref. 13) give algorithms for 
subdividing scalar displacements of one dimensional noise and 
two dimensional polygons. The subdivision of a fractal curve 
at its midpoint is also known as the midpoint displacement 
technique [ref. 1, pp. 233-234]. 

The actual first use of the midpoint displacement 
technique was for creating the Koch snowflake. The Koch 
snowflake uses an initiator and a generator [ref. 1, pp. 39- 


43] and replaces each initiator by the generator during each 


step in the recursive process. The displacement of the 
midpoint for the Koch snowflake was predetermined so that the 
generator had a specified shape. The method of displacing the 
midpoint is the most critical aspect of the midpoint 
displacement technique. There are three variations suggested 
for the midpoint displacement technique described in [ref. 1, 
pp. 233-234). They are the transversal midpoint displacement, 
the isotropic directions displacement and the random lengths 
displacement. The transversal midpoint displacement method 
allows the displacement of the midpoint to the left or to the 
right along the initiator according to specific rules as one 
moves away from the center. The isotropic displacement 
directions method randomizes the displacement directions of 
the midpoint. That is, if the line lies along the x axis, the 
displacement direction would be either positive y (up) or 
negative y (down) i.e., orthogonal to the x axis and to the 
Z axis. The random displacement lengths method allows the 
length or magnitude of the displacement to be random. Thus, 
the distance from the midpoint to the actual displaced point 
varies randomly. Using both the isotropic displacement 
directions method and the random displacement lengths method 
together affords complete randomness in the use of the 
midpoint displacement technique and is termed the random 
displacement method. 

The random displacement method, as described by Peitgen 


(ref. 14, pp. 133-136], is produced by the use of two 
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components. These are random variables and random functions. 
The random variables are used in the random functions and are 
created by a random number generator. The random function 
resulting from the random variables that are used, as Peter 
Sorensen [ref. 15, p. 160] indicates, achieve a displacement 
in the y direction. The y displacements are either up or down 
in the screen display, denoting positive or negative 
displacement, respectively, and are thus applicable in the 
creation of the surface of fractal mountains. When creating 
fractal mountains, the y direction displacements yield the 
altitudes of the various peaks and valleys. Care must be 
employed in the generation of the random numbers to prevent 
the random variables from becoming completely chaotic. 
Richard Voss [ref. 5] states that the random displacement 
method is a recursive generating technique first applied to 
the study of normal Brownian motion as early as the 1920’s by 
N.Wiener. In Voss’s explanation, the midpoint values for the 


random midpoint displacements are determined by the equations: 


Wi /2) 


0.5(V(0) + V(1)) + As (equation 2.1) 


Vv (72) 0.5(V(0) + V(-1)) + 42 (equation 2.2) 


where V(0) and V(1) and V(0) and V(-1) are the endpoints of 


the two initiators and 4, and 45 are gaussian random variables 


with zero mean and standard deviation 1. 
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Equations 2.1 and 2.2 yield the displaced midpoints of the 
two line segments V(0)V(1) and V(-1)V(0) as illustrated in 
Figure 2.1. The initial portion of the equations determines 


the midpoint of the two line segments. The 4's are the random 


adh 





Figure 2.1 Representation of Equations 2.1 and 2.2 


variables that determine the size of the displacement of the 
midpoint. 

Equations 2.1 and 2.2 are Similar to the equations used 
in the sample algorithms given by Fournier, Fussell and 
Carpenter [ref. 13] for the scalar displacement of one 
dimensional noise and of two dimensional polygons. 

The random midpoint displacement technique described above 
was first developed in one and two dimensions. The technique 
was expanded to three dimensions by Peitgen [ref. 14, pp. 95- 


105, 244-254]. Peitgen used the extension to three dimensions 
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to create three dimensional surfaces of fractal mountains. 
The fractal surface so created and displayed on a computer 
screen gives a realistic simulation of three-dimensional 
mountains as described by Patrick Bass [ref. 16]. A true 
three dimensional fractal mountain must have x, y and z values 
given for each location. A simple diagram of connections 
between the vertices of each small triangle by line segments 
from one base triangle to adjacent triangles can be easily 
programmed aS was done by Michiel van de Panne [ref. 17]. 
Alex Pentland [ref. 18, 19] interprets the third dimension as 
determining roughness or smoothness and isotropic or 
anisotropic conditions of the surface. The third dimension 
effect can also be described as modeling the texture of the 
fractal surface. Pentland (ref. 18] describes two ways that 
the surface texture is affected in image texture. These are 
1) projection foreshortening, a function of the angle between 
the viewer and the surface normal and 2) perspective texture 
gradient that is due to the increasing distance between the 


viewer and the surface. 


B. BASIC THEORY 

In this section the basic framework, the random midpoint 
displacement technique, is discussed. The termination 
criterion is discussed in the second section. The third 


section explains the three dimensional geometry and the last 
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section discusses the random variables used to displace the 


midpoints in the midpoint displacement technique. 


1. Midpoint Displacement Technique 
We return to the study of the recursive subdivision 
procedure. The random midpoint displacement technique forms 
the framework for the analysis. To review, the random 
midpoint displacement technique begins with a line segment and 
its two endpoints (also called vertices) oriented in three 
dimensions. The midpoint of the line segment is determined 


by adding the two endpoints for each of the x, y and z values 


separately and dividing by two. The second step is to 
randomly displace the midpoint in the y direction. The 
coordinate system shown in Figure 1.2 is used. The y 


coordinate is thus viewed as defining the height or vertical 
direction. The resulting new value for y with the previous 
unchanged values for x and z define the new midpoint. The new 
midpoint forms two new line segments with the two previous 
endpoints each of half the length of the previous line 
segment. These are each taken in turn as line segments for 
the next recursive call to the midpoint displacement routine. 

To develop a subdivision equation in three dimensions, 
one starts with a triangle. Three midpoints are determined 
with one on each side of the triangle. As seen in Figure 
2.2, one starts with a base triangle given by (P1,P2,P3), with 


each pair of vertices forming a line segment. From the two 
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endpoints for each edge, the midpoints M;,, M2 and Mz are found. 
By connecting lines between the midpoints M,, M> and Msz, 
smaller triangles are created, in this case triangles 1, 2, 
3 and 4. These new triangles are then recursively operated 
on separately by the midpoint displacement technique. Thus 
recursion is implemented on each of the triangles 1, 2, 3 and 
4 and the same subdivision technique is applied to each of 


the triangles in turn. 





Figure 2.2 Midpoint Displacement Technique 
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2. Determining the Termination Criterion 

The recursive process described in the previous 
section is continued until some termination criterion is 
satisfied. The termination criterion is at the control of the 
programmer and can depend on limitations imposed by the 
computer used in implementing the recursion. These 
limitations can also be subject to dynamic memory and disk 
space, collectively called memory, and screen resolution. 
Memory requirements imposed by the process can become quite 
large because each level of recursion creates four triangles 
from each triangle at the previous level. The more vertices 
that have to be stored further increase memory requirements. 
Since there is a limit to system memory, there is also a 
consequent limit on the number of recursive levels in the 
process. The amount of memory is known, so one can determine 
how many triangles and vertices can fit into the amount of 
space available by a combinatorial analysis. We can also 
count the numbers of triangles, edges and vertices developed 
in the recursive process. Each of these classes of objects 
play a part in determining how many iterations are possible. 
Each of the triangle, edge and vertex calculations has a 
separate defining recursion equation for the number of 
triangles, edges or vertices occurring. The equations are 
interrelated. For the following equations, the index i 


denotes the ith iteration. 
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For each new iteration, the number of triangles is 
four times the number of the previous iteration. Thus, the 
recursive equation for the number of triangles (T) created is: 

Tis) = 4 * T; i (equation 2.3) 
with the initial value Tọ= 1 

For the number of edges (E) created, the recursive 

equation is 
Eu = 2 * Ej + 3 * T; (equation 2.4) 
with the initial value Ej = 3. 

For each new iteration, we have double the previous 
edges plus three times the number of triangles of the previous 
level. This follows since, for every new recursion level, 
each edge becomes two edges and inside each triangle three new 
edges are generated when the midpoints are connected. There 
are no edges stored, per se, but the number of edges is needed 
to determine how many vertices will have to be stored for 
subsequent iterations of the process. 

The number of vertices (V) created is determined by 
the following recursive equation: 

Vier = Vi + E; (equation 2.5) 
with initial value Vp) = 3. 

The equation for the new number of vertices is 
dependent on the number of vertices of the previous iteration 
plus the number of edges for the previous iteration. On each 
of these edges, the midpoint process adds an additional vertex 


to the count at each recursive step. The resulting number of 
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vertices Vi; must be multiplied by three to determine the 
storage required since each vertex is described by a vector 
of three floating point numbers. Since each floating point 
number requires eight bytes of storage, the total number of 
bytes of storage needed is 24 times the total number Vi, of 
vertices. The number of bytes is used to determine the 
minimum amount of storage needed for each iteration or 
recursion level. This ultimately limits the number of 
recursive levels that can be processed. The maximum value V 
is determined by the minimum amount of available memory. 

To solve the above recursive equations, 2.3-2.5, the 
initial values for the first iteration used are T, = 4, E,; = 
9 and V, = 6. The process of determining how many triangles, 
edges and vertices are formed, can be determined at each 
iteration. The values for the first few iterations are given 
in Table 2.1. The closed form solution for the number of 
triangles is obtained from: 


To = 1 

T, = 4T, = 4 

T, = 4T, = 4%4 
It is not difficult to see that the general or closed form 
solution is given by 

T; = 4' (equation 2.6) 

By inspection of Table 2.1, one can see that 

E; = T; +V; -1 (equation 2.7) 
This relationship makes the development of the closed form 
solution easier. Substituting equation 2.7 into equation 2.5 
one has 

Vi = 2V; + T; - 1 
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Successive application of the recursion relationship yields 
Vie = 2(2Vj9 + Tiy -1) +7; - 1 
T FT ESEN OS 
Van (ee rg a tt OT. + Tf. ol - 2 
A ID o ET, MT - 1.- 2 - 4 
e AA A A A Dad 
ee + Te -1=224=8 


Evaluating this for i = 3 one obtains 
A O =a es 
The progression -1 -2 -4 -8 is a geometric 
progression and can be written as 1 - 2°, Leaping to a 


general formalation for V; we have 
i-t 
VS = 2'Vo + 2 2’ Tei-n-; Ze an 
J= 
By equation 2.6 we rewrite this as 


too, te 

Vi = 2'Vo +. Da En 
J= 

: 1-1 3 : u : 

Vi = 2'Vo + & Dar, 
J= 


Dr 2212 E + 1-2 
eo (2 = 1/21 0)) += 122 

ae M a y = 2! 

V; = aye u = o) al (equation 2.8) 
Substitution into equation 2.7 yields 

E; = a DIN, Oe an Byer) (equation 2.9) 

Equations 2.6, 2.8-2.9 gives the capability of 

determining how many triangles, vertices and edges, 
respectively, are generated at each level without needing the 


number at the previous level. The greatest number of levels 


of recursion that can be performed in the creation of fractal 
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mountains in this study is seven. With seven iterations the 
number of triangles is Ty = 16,384, the number of edges is E, 
= 24,768 and the number of vertices is V, = 8385. Since each 
vertex is determined by a vector of three components, x, y 
and z, and each component requires 8 bytes of storage. One 
has at the seventh level a memory requirement for vertices of 
24 X 8,385 or 201,140 bytes. Thus, slightly over .2 megabytes 
of storage is needed for storing the vertices of the base 


triangles that have been generated to the seventh level. 


Table 2.1 NUMBER OF TRIANGLES, EDGES AND VERTICES FOR 
RECURSIVE EQUATIONS TO I = 7 











Iteration Triangles Edges Vertices 

i De E V; 

0 1 3 3 
1 4 9 6 
2 16 30 15 
3 64 108 45 
4 256 408 1:53 
5 1024 1584 561 
6 4096 6240 2145 
7 16384 24768 8385 


The stopping rule for the recursive process can be 
terminated optimally by stopping at the iteration that results 
in the best use of the available memory. The rule is then 
applied to the equations 2.3-2.5. 

A second determining factor that could also affect the 
stopping criterion is a limitation on the size of the smallest 


addressable element on the output screen. The midpoint 
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displacement technique continually creates smaller and smaller 
triangles. The size of the triangles converge down to the 
smallest addressable quantity on the graphics screen, the 
pixel. The number of pixels on the screen will thus limit the 
total number of possible triangles appearing on a given 
computer display. If fewer iterations are performed, but the 
same size of the final triangle is desired, then the original 
base triangle must be set initially at a smaller size. The 
size of the triangle is measured by the number of pixels 
appearing on its edge. The larger the number of pixels 
available on a screen, the more triangles that can be created. 
Thus for a screen resolution of 1024 x 1024 pixels, it is 
possible to perform more levels of recursion then for a 
smaller screen resolution of say 640 x 320 pixels. 

Thus, for the ultimate stopping criterion, a 
combination of the above factors, available memory and screen 
resolution is needed. The greater the number of iterations 
performed by the midpoint displacement technique, the more 
natural the appearance of the fractal mountain. Therefore, 
one is motivated to use the computer in an optimal way, i.e., 


performing as many iterations as possible. 


3. Three Dimensional Geometry Based on Y Displacement 
To make a realistic fractal mountain, the mountain 
created must be viewed in three dimensions. The image in 


three dimensions is created by the random displacements of 
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the midpoint occurring along the normal to the x-z plane, i.e. 
in the y direction. Changing the values in the y direction 
creates the appearance of upward growth as would be desired 
for a fractal mountain. Since the midpoint displacement is 
done only normal to the x-z plane, fewer calculations are 
needed than would be necessary if x,z displacements were also 
included. A further reason for displacing only along the 
normal to the x-z plane is the appearance of the gapping 
problem which is discussed below in Chapters III and IV. The 
gapping problem refers to breaks in the polygonal covering 
that represents the topography of the fractal mountain. By 
changing the displacement only in the y direction, the 
resulting gapping problem, although still difficult and 
important, is minimized and becomes manageable. 

The displacement in the y direction helps create the 
natural looking fractal mountains surface. The appearance of 
the resulting surface is further determined by how rough or 
smooth a texture the surface has. In this study, a variable 
called ‘texture’ determines the smoothness of the mountain. 
This texture variable permits different mountain scenes, from 
rolling hills to rocky mountains as seen in nature, to be 


created on the computer screen. 
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4. How the Random Variables are Created for Displacing 
the Midpoints 


We wish the amount of displacement performed at each 
iteration of the midpoint algorithm to be determined at 
random. In this section, we iS the mathematical 
algorithm used to obtain the randomness in the midpoint 
displacement technique. 

The displacement of the midpoints in the y direction 
is based on random variables that are created with the help 
of a process of choosing random numbers drawn from a normal 
uniform distribution. The random number generator used 
produces a series of uniform random numbers that are then 
normalized to lie between zero and one. The normalized random 
numbers are drawn from a standard uniform distribution with 
mean O and standard deviation 1. An initial seed value is 
then chosen for the random number generator. The initial 
random number generator seed chosen for this study is an 
arbitrary value of 475836. 

The standard uniform random variables are converted 
to a Gaussian distribution and are placed into an array called 
RAND in the order that they are generated. The final array 
RAND is a table of 500 Gaussian distributed random numbers 
with mean zero and a standard deviation of one. The random 
number chosen for the displacement from array RAND is 
determined by a variable "loc". Loc determines the value for 


the number chosen from array RAND and ranges from 1 to 500. 


23 


The value of loc is important in the algorithm due to 
the two fundamental problems that occur in achieving natural 
looking fractal mountains. These two fundamental problens of 
granularity and gapping are described and defined in the next 
chapter. Following that description, in Chapter IV a 
description of the impact of the granularity and the gapping 


problems on the algorithms used in the implementation of the 


random midpoint displacement technique is made. 
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III. PROBLEMS AFFECTING REALISM 


Two problems affect the realism we see in the creation of 


fractal mountains. The problens are granularity and gapping. 


A. GRANULARITY 

A major concern in creating fractal mountains is the 
realism of the resulting picture of the mountain. The 
appearance of a natural look in the picture of a mountain is 
determined in part by the texture given by the algorithm that 
creates the mountain. The extremely regular pattern 
consisting of rows of valleys and ridges shown in the fractal 
mountain in Figure 3.1 is not representative of actual 
mountains. This regular pattern is called granularity. The 
more patterned or granular a fractal mountain is, the less 
realistic the mountain appears. To eliminate the granularity 
as shown in the fractal mountain of Figure 3.1, more 
randomness must be used in the midpoint displacement 
technique. The more random the midpoint displacements 
produced by the algorithm, the less granular will be the 
appearance of the fractal mountains. For example, the 
technique used to produce the mountain in Figure 3.1 used only 
8 different random variables. The methods used to minimize 
the granularity by the introduction of more randomness are 
discussed in Chapter IV. Additionally, the 
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granularity is affected by the orientation of the triangles 
An important factor affecting the granularity is how the 
subroutine calls are made for each triangle being operated on 
by the midpoint displacement technique. 
B. GAPPING 

Another problem that affects the natural look of the 
fractal mountain is the "gapping problem". The gapping 
problem occurs when one applies the midpoint displacement 
technique to the common side of two adjacent triangles. 


Referring to Figure 3.2a, the two adjacent triangles are A and 





Figure 3.1 Example of Granularity 





B. They share the common side P,P, with the midpoint indicated 
by M. The midpoint displacement technique displaces M first 
in triangle A and then in triangle B. In Figure 3.2b the point 
M has been displaced in triangle A to point M, and in triangle 
B to point M. The resulting polygon P3M,P,M, defines a gap or 
discontinuity in the figure. In order to eliminate the gap, 
the vector displacement of M in triangle A must be the same 
as the vector displacement of M in triangle B. Chapter IV 


discusses the technique used to insure displacements without 


a gap. 





Figure 3.2 Gapping Problem 
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C. INTERRELATIONSHIP BETWEEN GRANULARITY AND GAPPING 

In order to eliminate the problem of granularity in the 
creation of fractal mountains, randomness is needed in the 
midpoint displacements, which also results in a realistic 
mountain. On the other hand, the problem of gapping from two 
different midpoint displacements of the common side of 
adjacent triangles demands that the two displacements of a 
midpoint on a common side be equal. This means a certain 
regularity has to be imposed on the midpoint displacement 
technique. A conflict results between the solution of the 
granularity problem and the solution of the gapping problen. 
Thus, a trade off has to be made between randomness and 
regularity . The effect of granularity and that of gapping and 
the combination of the two problems is further discussed in 
Chapter IV. The simultaneous solution, finally adopted for 
this study, for both the granularity problem and the gapping 


problem is presented in Chapter V. 
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IV. PUTATIVE ALGORITHMS FOR THE GRANULARITY AND GAPPING 


PROBLEMS 


In Chapter III, the problems of granularity and gapping 
affecting the naturalistic look of a fractal mountain were 
defined and briefly discussed. This chapter discusses the 
algorithms that have been considered to solve the granularity 
problem and the gapping problem. The various approaches 
discussed here did not lead to a simultaneous solution to both 
of the problems. The discussion of these algorithms is 
included because it seems useful to do so. The discussion 
shows the complications that can occur in an effort to create 
natural looking mountains. The solution that simultaneously 
handles the granularity problem and the gapping problen, 
incorporating some of the approaches used in this Chapter, is 


presented in Chapter V. 


A. RANDOM VARIABLE EFFECTS ON GRANULARITY 

The granular appearance of the fractal mountain depends 
on the randomness of the midpoint displacements. It is 
possible the resulting granularity is caused by some 
undiscovered regularity in the random variables generated by 
the random number generators. To determine if the granularity 


is due to the quality of the random variables, the Gaussian 
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distribution random number generator, described in chapter II, 
was replaced by a linear congruential random number generator. 
The linear congruential method creates a sequence, ati), 

of random numbers by using the following equation: 
a(i) = ((a(i-1) * b) + 1) mod m (equation 4.1) 
An arbitrary initial value a(0) is needed to start equation 
4.1. Constants b and m for the multiplier and the modulus, 
respectively, are also used in the process. Each different 
initial value yields a different random sequence. The next 
random number in the sequence is obtained from a(0) by 
multiplying a(0) by the constant multiplier b and adding 1. 
The operation [(a(0)*b)+1] is performed modulo the second 
constant m. The resulting value is taken as the second random 
number a(1). The process continues, producing 
a(2)., atinado The constants m and b are chosen 
subject to some constraints. Sorensen (ref. 20] states that 
the best choice for m is a power of 10 or a power of 2 
depending upon the index of the number system used. Further 
he states that b should be a number one digit less than m. 
An additional constraint for b is that it be of the form x21 
where x is an even number. For example with m=1024, and with 
X arbitrarily chosen as 8, b becomes 821. (notice m has four 
digits and b has three digits). The resulting random number 


sequence is composed of integers between 0 and m-1. For this 
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random number sequence to be used in the midpoint displacement 
technique, the numbers need to be normalized to lie between 
(oij: l 

Using the random numbers produced by the linear 
congruential method instead of those produced by the normal 
uniform distribution method, discussed in Chapter II, does not 
materially change the granular texture observed in Figure 3.1. 
It is therefore concluded that the granularity problem is not 
due to the quality of the numbers produced by the random 
number generator. The normal uniform distribution generator 
is therefore retained in the actual implementation. It is 
therefore further concluded that the granularity problem is 
the result of some different aspect of the implementation of 
the midpoint displacement technique. 

To understand the effect of the random variables on 
granularity, the effect of the selection of the variable loc 
was investigated. The variable loc is used to choose the 
random number chosen from the array RAND. To begin the 
process, loc can be arbitrarily chosen equal to 1. The first 
random variable used is thus the random variable in the first 
position of array RAND. When the base triangle P1 P2 P3 is 
operated on by the random midpoint displacement technique, the 
midpoints Pmidl, Pmid2 and Pmid3 are created (see Fig 4.1). 


The first midpoint displaced is Pmidl, so the displacement of 


S1 


Pmidi in the y direction uses only the first number from RAND 
(iver, loc = 1). A predetermined value of the variable 
texture, as described in Chapter II, is also incorporated in 
the equation: 
Pmidl(y) = (texture * RAND(loc)) + Pmidl(y) (equation 4.2) 
For each of the succeeding sides in the triangle, the 
pointer loc is successively incremented by 1. One obtains 


the equations: 


Pmid2 (y) = (texture * RAND(loc+1))+Pmid2(y) (equation 4.3) 

Pmid3(y) = (texture * RAND(loc+2))+Pmid3(y) (equation 4.4) 
Note: the hatted variables are updated versions of the 
previous variables. The hats shall be suppressed in the 
sequel. 

The random displacements of the midpoints is a recursive 
process. We also update the random numbers drawn from RAND. 
At the next level of recursion, we increment the value of the 
pointer in the table RAND (i.e., loc = loc + 1). For the next 
set of midpoint displacements, the pointer loc is set equal 
to the second, third and fourth positions of RAND. This 
process is repeated through each level of iteration (See Fig 
As) The positions from RAND used for each midpoint are 
denoted by arabic numbers. For six iterations, there are up 
to eight different random variables chosen from array RAND. 


In this case, the random variables are those in locations 1 
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to 8 of the array RAND. The process shown which makes smaller 
triangles is also performed for triangles II, III and IV. 
Other smaller triangles are created ont i the triangles are 
as small as a predetermined size. This yields a structure 
such as seen in the top triangle of Figure 4.1. The regular 
pattern of the numbers as seen in Figure 4.1, i.e., the 
selection of locations of RAND indicates that the midpoint 
displacements are not truly randon. This method of selecting 
loc leads directly to the granularity of Figure 3.1. Thus, 
somehow the regularity in the selection of the random 


variables needs to be avoided. 
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Figure 4.1 Midpoint Displacement Technique's Seed Values 
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B. MODIFICATION OF THE MIDPOINT DISPLACEMENT TECHNIQUE TO 
ELIMINATE GRANULARITY 
As the random number generator has been shown to not have 
a deleterious effect on the granularity problem, other 
solutions for eliminating the granularity need to be 
considered. Both the effect of weighted displacements of the 
midpoint and the effect of small based triangles were 
considered in order to eliminate granularity. 
1. Weighted Displacements of the Midpoints 
The first modification to the midpoint displacement 
technique considered is to create weighted displacement of the 
midpoints. The term weighted displacement is used to describe 
different maximum displacements of the midpoints as a function 
of the size of the initiator (the line segment) at each step. 
When the size of the initiator is smaller than a predetermined 
value, the maximum displacement of the midpoint is reduced in 
the weighted midpoint displacement technique. The original 
implementation of the midpoint displacement technique used in 
this study, was taken from work of Michael Gaddis [ref. 21] 
and is described above as the original implementation. This 
implementation of the midpoint displacement technique resulted 
in the granularity shown in Figure 3.1. A detailed 
description of this algorithm is needed in order to understand 


fully the modification of weighted midpoint displacements 
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developed for solving the granularity problen. 
Original Algorithm (Equations 4.5-4.10) 

To determine the displaced midpoint in the first 
implementation of the midpoint displacement technique, the 
midpoint of the y values between two endpoints of the 
initiator is first found. 

Initialization Equations (Equations 4.5-4.7) 

The equations to determine the midpoints for the three 

edges of the triangle first find the midpoint of the x, y and 


z values from the endpoints. The equations are: 


Pmidl(x) = (P1(x) + P2(x))/2.0 (equation 4.5) 
Pmid2(x) = (P2(x) + P3(x))/2.0 (equation 4.6) 
Pmid3(x) = (P3(x) + P1(x))/2.0 (equation 4.7) 


Similar equations are usedto find the three midpoints 
in terms of their y and z coordinates. Once the midpoint on 
each edge is found, the y coordinates of the midpoints are 
randomly displaced. 

Displaced Equations (Equations 4.8-4.10) 
The random displacement for the three midpoints in the 


y direction are given as: 


Pmidl(y) = (texture * RAND(loc)) + Pmidl(y) (equation 4.8) 
Pmid2(y) = (texture * RAND(loct1))+ Pmid2(y) (equation 4.9) 
Pmid3(y) = (texture * RAND(loc+2))+ Pmid3(y) (equation 4.10) 


Texture (described in Chapter II.B.3) is the roughness 


factor and takes on a value in a range between 1 and 20 to 
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give the roughness or smoothness of the surface. RAND(loc) 
is explained in Chapter II.B.4 and is the random number used 
to determine the amount of displacement of the midpoint value. 
Note three different locations are chosen for the three 
equations. The displacement of the midpoint of each edge, 
Pmid1(y), Pmid2(y) and Pmid3(y) are then determined by the 
equations 4.8, 4.9 and 4.10, respectively. 
Weighted Displacement Algorithm (Equations 4.11-4.16) 
In the modified algorithm, where weighted displacement 
of the midpoints is attempted, the initial equations 4.5, 4.6 
and 4.7 are retained, but the displacement equations 4.8, 4.9 
and 4.10 are altered. The difference in the equations used 
is determined by the distance between the two endpoints of the 
two points in the x and z coordinates, measured in pixels. 
For each of the endpoints, the distance is determined. The 
equations are: 
Selection Equations (Equations 4.11-4.13) 
Pabl= sgqrt[ (P1[x]-P2[(x])* + (Pl[z]-P2[z])* ] (equation 4.11) 
Pab2= sqrt[ (P2[x]-P3[x])* + (P2[z]-P3[z])* ] (equation 4.12) 
Pab3= sgqrt((P3[x]-P1[x])% + (P3[(2]-P1(2])% ] (equation 4.13) 
Pabl, Pab2 and Pab3 are used in the criteria to 
determine whether the displaced equations are to be modified. 
If any of the values Pabl, Pab2 or Pab3 are greater than a 


selected value €, then the respective equations used are the 
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original displaced equations 4.8, 4.9 and 4.10 for the given 
value. The use of equations 4.8-4.10 allow a large 
displacement in the y direction, in general. As the distance 
between the points measured in the x and z coordinates gets 
smaller, the possible displacement needs to be smaller than 
when the size of the initiator is large. This change in 
displacement prevents the possibility of a valley and a peak 
in a mountain scene being right next to each other. This can 
scour naturally, fer example, in a* rift valley setting. Por 
now we seek more ordinary mountain settings. 

If, on the other hand, any of Pabl, Pab2 or Pab3 are 
less than or equal to e, then the respective equations used 
for displacement are a function of the distance of the x and 


z coordinates: 


Modified Displaced Equations (Equations 4.14-4.16) 
Pmidl(y)=(RAND(loc) *texture* (Pabl/e)) + Pmidl[y] 
(equation 4.14) 
Pmid2(y)=(RAND(loct1) *texture* (Pab2/e)) + Pmid2[y] 
(equatuon 4.15) 
Pmid3 (y) =(RAND(loc+2)*texture* (Pab3/€)) + Pmid3[y] 
(equation 4.16) 


Mete, e is set to 50.0 and loc is the location of the 
random number chosen from the array RAND. 
e is set to 50.0 to allow the number of edges of each 


triangle to occur at different recursion levels if the 
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triangle is an acute, an equilateral or an obtuse triangle. 
The ending criteria (see Chapter VI.D) is set at either 16 
pixels or 9 pixels for the length of an edge of a triangle. 
Because edge lengths are successively halved at each 
iteration, for the edge of a triangle to be less than 16 
pixels, the resulting edge length would be strictly between 
8 and lee Suppose the edge length when the stopping criterion 
is met occurs at the fifth recursion level. Then at the 
fourth recursion level there would have to be an edge of 
length between 16 and 32. The third recursion level would 
have edge lengths between 32 and 64. The value of 50.0 lies 
approximately 2/3 of the way from the value of 32 to the value 
of 64. Therefore, if the base triangle is chosen to have 
edges of length between 250-500 pixels, as is common in this 
study, then the following will occur for the three types of 
triangles. If the triangle is an acute triangle when the 
third recursion level is completed, one of the three edges 
would be under 50 pixels in length and therefore the weighted 
displacement equation would be used in the next iteration on 
that edge. Then on the fourth recursion level the other two 
edges would also use the weighted displacement equations. If 
the triangle is an equilateral triangle, the three edges of 
the triangle would all use the weighted displacement equations 


at the same recursion level. If all sides are less than 350 
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pixels in length, the recursion level at which the weighted 
displacement equations are used would be the third level. If 
the triangle is an obtuse ICA then at the third 
recursion level two of the three edges would use the weighted 
displacement equations. The third edge would use the weighted 
displacement equation at the fourth recursion level. The 
results for the acute, equilateral and obtuse triangles makes 
e equal to 50.0 an attractive choice. 

If, for example, only Pab1 is less than e, then 
equation 4.8 is replaced by equation 4.14 and the equations 
4.9 and 4.10 are used in the other two respective cases. 
Corresponding changes for Pab2 and Pab3 are made if these 
numbers are less than e. 

Using the weighted displacement algorithm for the 
midpoints gives only a slight improvement in the realism of 
the fractal mountain. Although there was less granularity in 
the picture corresponding to Figure 3.1, the fractal mountain 
still did not look very realistic. The weighted displacement 
technique of the midpoints is thus not the final solution for 
improving the realism of the fractal mountains, so further 
modifications were sought. However, since there was some 
improvement with the weighted displacement technique, this 
technique is incorporated in the implementation used to 


achieve the naturalistic effect in fractal mountains. 
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2. Smaller Base Triangles 

The next attempt to remove the granularity is to use 
smaller base triangles. The original implementation began with 
a base triangle with a couple hundred pixels in length for 
each edge the size seen in Figure 3.1. The original base 
triangle is then divided into four arbitrary sized triangles 
of about 40-50 pixels for each edge as seen in Figure 4.2. 
These four new triangles are then used as base triangles in 
the algorithm implementing the midpoint displacement 
technique. The resulting fractal mountain using the smaller 
base triangles is seen in Figure 4.3. The four triangles each 
have a different granular appearance. On the right, where 
triangles III and IV join, small gaps can be seen. There are 
also gaps between each of the pairs of triangles, but on only 
the one common side of adjacent triangles are the gaps large 
enough to be seen without magnification. Each triangle still 
exhibits a distinctive granular appearance. Because there are 
four different granularities, the fractal mountain looks more 
realistic compared to the original fractal mountain as shown 
in Figure 3.1. This solution of the granularity problem with 


small base triangles, however, introduces a gapping problem. 
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C. WHY THE ORIGINAL IMPLEMENTATION DID NOT RESULT IN GAPS 

In sections 4A and 4B, the ar different random 
variables and smaller base triangle sizes were considered to 
deal with the granularity problen. In section 4B, smaller 
triangles helped solve the granularity, but introduced gaps 
into the solution of the midpoint displacement technique. In 
this section, the gapping problem is further considered. The 
fractal mountain of the original implementation seen in Figure 
3.1 contains no gaps. Although the gapping problem did not 
exist within the original implementation of the midpoint 
displacement technique, we have seen that efforts to create 
better granularity can introduce a gapping problem. A 
description of the algorithm in the original implementation, 
may help in an understanding of the modifications needed for 
the algorithm to solve both the granularity problem and the 
gapping problem simultaneously. 

As discussed in Chapter II, the creation of fractal 
mountains is accomplished recursively. The recursion process 
is fundamental for the creation of fractal mountains. But 
just as important, for realistic fractal mountains is the 
implementation of the recursion process. In the original 
implementation, the subroutine used to implement the midpoint 


displacement technique is called mountain generate. By the 
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Figure 4.2 The Four Base Triangles 








Figure 4.3 Fractal Mountain using Small Base Triangles 
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process of recursion, the mountain_generate algorithm is 
called from within mountain_generate. As seen in Figure 4.4, 
Triangle P1 P2 P3, when aaa by the subroutine 
mountain_generate, is divided into four new smaller triangles 
within the original triangle. Each of the new triangles forns 
the basis for a new level of recursion. The process continues 
until some smallest acceptable triangle size is reached by the 
stopping criterion. This results in the triangle structure 
such as seen in the top triangle of Figure 4.1. In Figure 
4.4, the new midpoints Pmidl, Pmid2 and Pmid3 when connected 
create the four new triangles I, II, III and IV. Notice 
triangle I is in the upper part of the original base triangle 
P1 P2 P3. The orientation of each triangle I, II, III and IV 
is important in achieving the desired end of eliminating the 
gapping problem. The triangles have to be rotated in an 
appropriate way for each call to mountain generate. 

To implement the recursive call for the midpoint 
displacement technique on each of the triangles, I, II, III 
and IV, a call to mountain_generate is used. The 
mountain_generate call uses the three endpoints of the 
respective triangle as input variables and a fourth input 
variable, for indicating the location in RAND for the value 
for the random number used. The order of the three vertices 


of the triangles in the call is crucial. Referring to Figure 
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Figure 4.4 Seed Values Used for the Recursive Process 


4.4 for each triangle, the following recursive Calls to 
mountain generate are used. 
Recursive Calls 
mountain _generate(P1,P2,P3,10c); for Base Triangle 
loc =#locet ai; 
mountain generate(Pmid1,P2,Pmid2,loc); for Triangle I 
mountain generate(Pmid3,Pmid2,P3,loc); for Triangle II 
mountain _generate(P1,Pmid1,Pmid3,loc); for Triangle III 
mountain generate (Pmid2,Pmid3,Pmidl,loc)'; 
for Triangle IV 
As can be seen from Figure 4.4, the above recursive call 
for triangle I starts with the left endpoint Pmidl as the 


first variable, the uppermost endpoint P2 for the second 


variable and the right most endpoint Pmid2 for the third 
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variable. The first, second and third endpoints are 
determined as the corresponding first, second and third 
variables sent to the routine a ene sues The order 
of the input variables to mountain generate is important, as 
this order determines which value of loc is used on each edge 
of the triangle for displacing the midpoints. The second, 
third and fourth recursive calls rotate the triangles II, III 
and IV so that appropriate loc values are used on the common 
edges of adjacent triangles. 

Rotation of the triangle is accomplished by the order in 
which the endpoints are presented in the mountain_generate 
call. As indicated previously, triangle Pl P2 P3 forms the 
first level of recursion, while triangles I, II, III and IV 
are in the second level of recursion. For each level of the 
recursion, the same three consecutive values are used for loc. 
Thus, values of loc equal to 1, 2 and 3 are used in the first 
level of recursion. Loc is then incremented by 1. For every 
triangle created in the second level of recursion the values 
of loc are taken to be 2,3 and 4. These values are denoted 
by the numbers on the respective edges in Figure 4.4. 
Specifically, in the call for triangle I in Figure 4.4, loc 
will be 2 for the midpoint between the first endpoint Pmidl 
and the second endpoint P2. Similarly, loc is 3 for the 


midpoint between the second endpoint P2 and the third endpoint 
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Pmid2. The last edge of the triangle has loc equal to 4 used 
for displacing the midpoint between the third endpoint Pmid2 
and the first endpoint Pmidl. As can be seen in the recursive 
calls, the same values of loc i.e., 2, 3 and 4, are used three 
more times for each of the triangles II, III and IV in Figure 
4.4. The value of loc used for triangles II, III and IV is 
determined in a similar manner as that described for triangle 
I using the last three recursive calls shown above. 
Triangles I, II and III are oriented so that the common sides 
of the adjacent triangles have the same loc and therefore, the 
same displacement. This matching of common sides of adjacent 
triangles yields a solution to the gapping problem because the 
same random number is used in the displacement process for 
both calls of the algorithm on a common side. 

Notice however that only four different loc values are 
used. And each value is repeated four times at the given 
level of the recursion in the algorithm. Repeated use of the 
Same value of the loc implies that the displacement of the 
midpoints is not completely random as discussed previously. 
In addition, there is regularity or granularity introduced in 
the fractal mountain by the orientation pattern of the 
triangle. Although the original implementation has no gaps, 


it does tend to perpetuate the granularity problem. 
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In order to solve simultaneously the granularity problem 
and the gapping problem, careful attention to detail is 
required. Thus not only is the basic technique used to create 
fractal mountains important, but also the actual 
implementation of the midpoint displacement technique is 
important. The implementation is important in order to insure 
true randomness with no gaps in the creation of the fractal 
mountain. In the next chapter, the algorithms finally 
developed to solve both problems are given. That is, the gaps 


are controlled and an amount of granularity that can be called 


acceptable is achieved. 
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V. IMPLEMENTATION USED TO ACHIEVE A NATURALISTIC EFFECT IN 
THE CREATION OF FRACTAL MOUNTAINS 

This chapter explains the implementation that ensures a 
naturalistic look to fractal mountains. The first section 
explains how the random variables are used to achieve a 
naturalistic look. In section B, the hidden surface removal 
and the lighting calculations are discussed. In the last 
section, the algorithms dealing with the intersections of the 


outline for the mountain range are discussed. 


A. HOW RANDOM VARIABLES ARE USED TO CREATE A NATURALISTIC 

LOOKING FRACTAL MOUNTAIN 

Random variables are used to create in the fractal 
mountain the realism seen in nature. As was seen in Chapter 
IV, the modifications to the algorithms using the midpoint 
displacement technique showed that the technique was important 
in producing realistic mountains. Also, unless great care is 
taken with the ¡implementation of the technique, the 
simultaneous solution of the granularity and gapping problens 
cannot be achieved. As was discussed in Chapters III and IV, 
there is a trade off between the degree of randomness employed 
and the degree of regularity used. On the one hand, to solve 


the granularity problem, a great deal of randomness is needed 
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in the midpoint displacements. On the other hand, the gapping 
problem requires a certain amount of regularity in displacing 
the midpoints. 

The random variables used in the program are formed into 
a uniform normal distribution. To solve the gapping problen, 
the midpoint displacement on the occasions when the common 
side of two adjacent triangles are considered, must be the 
same. The two endpoints of the common side of adjacent 
triangles are the same and the values of x and z coordinates 
are not displaced by the algorithm. Therefore, the value of 
loc, used to determine the random variable to displace the 
midpoint in the y direction, is chosen to be a function only 
of the x and z coordinates of the two endpoints. In this way, 
the displacements of the same point in the two triangles is 
never different. Thus, no gaps can develop. 

The following equations (5.1-5.3) are used to determine 
the value of the location in RAND for each of the midpoints 
of the edges of the triangle. To determine the value of loc, 
the two x values and the two z values are added. Addition is 
used for the two endpoints instead of a subtraction so that 
the value is not dependent on the order in which the vertices 
are chosen. Since the value should range between O and 499 


the addition of the x and the z coordinates are modulo 500. 
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For displacing the midpoint Pmidl, as seen in Figure 4.4 


the equation is: 


PT = (P1[x]+P2[x])+(Pi[z]+P2[z]) 
loc = PT mod 500 (equation 5.1) 





Figure 5.1 Seed Values for the Recursion Process 


For displacing the midpoint Pmid2 the equation is: 


PT = (P2[x]+P3[x])+(P2[z]+P3[z]) 
loc = PT mod 500 (equation 5.2) 


and the equation for midpoint Pmid3 is 


PT = (P1[x]+P3[x])+(Pı[z]+P3[z]) 
loc = PT mod 500 (equation 5.3) 


Since there are 500 values in array RAND, the value of 
loc, must be between O and 499. Using the modular 


computation routine gives a good distribution of numbers from 


O to 499. 
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In Chapter IV, it was shown that in the original 
implementation the triangles in the iterative procedure were 
oriented a certain way to prevent the appearance of gaps. 
Equations 5.1-5.3 solve the gapping problem without having to 
turn the triangles for each level in any particular 
orientation as was needed in the original implementation. The 
mountain_generate calls for the recursion process are the 
following: 

mountain _generate(Pmid1,P2,Pmid2,ilev) 

mountain generate(P3,Pmid3,Pmid2,ilev) 

mountain generate (Pmid1,Pmid3,P1,ilev) 

mountain generate (Pmid1,Pmid3,Pmid2,ilev) 

Notice in the above calls to mountain _ generate that the 


triangles are presented to mountain generate in one order. 


Since the order does not matter, the endpoints above are only 


one order that can be used. All combinations of the three 
endpoints are possible. The variable ilev is used to 
determine the recursion level being executed. Ilev is used 


so that if there is more than one base triangle, all the base 
triangles will have the same number of recursion levels and 
the final iteration will result in triangles of the same 
approximate size. 

Although there is regularity in displacing the midpoint 
of the common side of two adjacent triangles, there is also 
enough randomness to prevent the patterns characteristic of 
the granularity problen. That is, the variable loc is 


computed in such a way to allow any of the values of RAND to 
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be accessed. The original implementation used only 8 random 
variables in processing through five levels of recursion. The 
present implementation uses all of the 500 random variables 
in RAND at least in theory. In the case where the recursion 
level results in more than 500 midpoint displacements, there 
will necessarily be repeated random variables. This 
repetition will not occur in any regular pattern and therefore 
will not result in regularity or granularity appearing in the 
fractal mountain. 

Another idea providing randomness is that of weighted 
displacement (Chapter IV.B). The idea of weighted 
displacement has been implemented to make the mountains as 
naturalistic looking as possible. 

With the simultaneous solution of the granularity and 
gapping problems to give the desired realism to the fractal 
mountains, the subject of hidden surface removal and of the 
mathematical calculations for lighting can be explained. 
These two subjects are important components in producing a 


natural look in the display of fractal mountains. 


B. HIDDEN SURFACE ELIMINATION AND LIGHTING CALCULATION 
Hidden surface elimination and sophisticated lighting 
routines can be used to give fractal mountains a natural look. 
Z-buffering is used to achieve the hidden surface elimination 
and the lighting is used to provide the effect of different 


times of day. 
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1. How Z-buffering Adds to Fractal Mountains 

Since fractal mountains are being created in three- 
dimensions in our study, there is a need to develop the 
concept of hidden surface elimination in the algorithm for the 
computer being used. The z-axis is portrayed as being 
perpendicular to the screen or to the xy plane. The concept 
used is called Z-buffering. [ref 22] Z-buffering adds to the 
naturalistic look by achieving the hidden surface elimination 
in the fractal mountain. 

2. Mathematical Calculations for Lighting 

Lighting and shading also affect the naturalistic look 
of fractal mountains. We use the IRIS workstation lighting 
capabilities. The creation of a fractal mountain requires the 
operator to provide a definition of 1) the color and the 
surface properties of the fractal mountain, 2) the color, 
location and direction of the light source and 3) the position 
and view direction of the observer. The surface properties 
for the fractal mountains are characterized by three factors. 
These are ambient light, diffuse light and specular light. 
Ambient light is a light of uniform brightness caused by 
multiple reflections from many surfaces of the fractal 
mountain and is non-directional. Light from ambient and point 
sources, which is scattered in all directions by the surfaces, 
is called diffuse light. It results from the surface 


roughness or graininess. Point sources create highlights 
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called specular reflection on shiny surfaces. Since specular 
light is mainly from smooth surfaces very little specular 
light is used for general fractal mountains. We note an 
exception in the next chapter. The ambient light and diffuse 
light vary according to the type of color and surface desired 
for the fractal mountain texture. Each of the three light 
types is specified by a red, a green and a blue component. 
Each component is given a value between O and 1. A forested 
mountain, for example, would have a large green component and 
small red and blue components. 

Characterization of the light source is an important 
factor in the lighting. The light source, which is taken 
to be the sun, consists of an ambient light component, a 
color component and the position of the light. The ambient 
light and color are specified by red, green and blue values 
between 0.0 and 1.0 as in the specification of the mountain 
surface light properties. The color of the sun is again 
specified by its red, green and blue components. The position 
of the light source or the sun is determined by a four 
component vector with values given by x, y, z and w. The w 
value determines if the light source is local or at infinity. 
If the light source is local, w is set equal to 1.0 and if w 
is set equal to 0.0, the light source is taken to be at 
infinity. When the light source is taken to be the sun, it 
is always placed at infinity. The x, y and z values 


determine the position of the light in the local case or the 
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direction of the light E when the light source is set at 
infinity. The x, y and z values for the location of the light 
source can be varied to achieve morning, noon or evening 
lighting effects. 

An additional property of the lighting consists of the 
overall level of light falling on the entire scene. This is 
referred to as scene ambient light. The scene ambient light 
parameter controls the general level of light in the fractal 
mountain scene. When the scene ambient light is low such as 
at morning, the whole scene is darker than when the scene 
ambient light is high such as at noon. To define ambient 
light, the three values of the red, green and blue components 
between 0.0 and 1.0 must be given. The location of the viewer 
is also a variable in the light calculations. To determine 
if the viewer is a local viewer or a viewer seated at 
infinity, a logical variable must be set for the program. 
This logical variable is 1.0 for TRUE if the viewer is local 
and 0.0 for FALSE if the viewer is at infinity. 

Once the above properties of the mountain surface, the 
light source and the scene light have been defined, the 
lighting calculations for the fractal mountain can be made. 
To make the calculations, the light perceived by the viewer 
depends on the angular relationships between the light source 
direction, the view direction and the surface normals. The 
light source direction is the position of the sun as specified 


in the program with respect to the surface of the mountains. 
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The view direction is the direction of the observer. The 
surface normal is the orientation of the surface of each 
triangle that is used in creating the fractal mountain. The 
normal is perpendicular to the surface of each triangle. 

This surface normal is a vector and is calculated from the 
vertices of the triangle by using the dot product, as shown 


in the following equations from the progran. 


NLS[0]=((P2[y]-P1[y])*(P3[2]-P1[2]))-((P3[y]-P1[y])*(P2[2]- 


P1[Z])) 
NLS [(1J)=((P2[Z)-P1[2]) *(P3[xX]-P1[x]))-((P3[2)]-P1[2]) *(P2[x]- 
P1[x])) 
NLS[2]=((P2[x]-P1[x])*(P3[y]-Plly]))-((P3[x]-P1i[x]) *(P2[y]- 
P1[y])) 


The surface normal vectors needs to be of unit length. 
Thus, x**2 + y**2 + z**2 equals 1. With x = NLS[O], y = 


NLS[1]) and z = NLS[2) the normalized surface vector is given 


by 

NL1[0]) = NLS[0]/sqrt((NLS[0]**2) + (NLS[1])**2) + (NLS[2]**2)) 
NL1(1)] = NLS[1]/sqrt((NLS[0]**2) + (NLS[1]**2) + (NLS[2]**2)) 
NL1[2] = NLS[2]/sqrt((NLS[0]**2) + (NLS[1}]**2) + (NLS[2]**2) ) 


When the viewer is positioned at infinity and the light 
source is at infinity, as is the typical case for fractal 
mountains, only one normal vector for each triangle needs to 
be given. The calculations for lighting the fractal mountain 
are all done automatically by the software provided with the 
IRIS workstation. Once the surface properties, the light 
source and the light model have been defined and the normal 


is sent to the workstation, the calculations for the color of 
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the triangle can be determined. When the light source 
position changes such as from morning to evening, different 
effects will appear on the fractal mountain which helps create 


a realistic image. 


C. MOUNTAIN RANGE OUTLINE 

The mountain range outline can be drawn anywhere within 
the window. Since one can draw anywhere within the window, 
there must be some routine that checks for intersections and 
duplicated y values for each x value. A description of this 
routine follows. 

Each point of the outline is connected by a straight line 
with the previous point. Each pixel with an x coordinate 
must have a unique y coordinate. To accomplish one y 
coordinate for each x coordinate, an interpolation must occur 
between pairs of points. The interpolation is done by first 
determining the slope between the two points. Starting with 


point (x,,y,) and point (x,,y,), the slope is computed as: 


sl = (Y, - y,) / (X> - Xy) (equation 5.4) 

Once the slope has been determined, the equations used to 
determine the y coordinate for each x coordinate are 
determined by: 

x; = i 


Y, = (sl*x,) + y, (equation 5.5) 


Here, i takes on values in the range of 0 to 1049. 
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Once the array y, is filled, the values computed and 
placed in the array are compared to all the values given when 
creating the outline. If there are duplicate y coordinates 
for a given x coordinate, then the larger value of y is used. 
The result of this operation is to remove any back tracking 
or any intersections. When one of these pairs of duplicated 
y values does occur, the largest value is taken for the 
outline of the mountain range. 

The algorithms for creating a realistic looking fractal 
mountain for use within the ISFM system have now been covered. 
The next chapter explains how these algorithms are actually 


used with the interactive system. 


58 


VI. USE AND IMPLEMENTATION OF THE INTERACTIVE SYSTEM 
FOR FRACTAL MOUNTAINS 

The Interactive System for Fractal Mountains (ISFM) is an 
interactive system to facilitate the creation of fractal 
mountains for an observer. ISFM creates the fractal mountains 
by means of an interface using pop-up menus. Pop-up menus 
permit the selection of options to change the type of 
mountain, to change the location of the light source and to 
create single mountains or mountain ranges. The interactive 
features permit the observer to see the mountains at different 
times during the day. 

Since the interactive system is controlled by pop-up 
menus, the menus give the different options that are 
available. The main menu is: 

Fractal Mountains 
Tutorial 

Texture 

Lighting Control 
Structural Factor 
Create a Mountain 
Outline 

Create a Mountain Range 
Clear and Reset 

Exit 

Since the system was done to give a systematic approach 
to the creation of mountains, two types of fractal mountain 


output have been developed. The first type, the single 


mountain option, 'Create a Mountain', has been included to 
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give a quick method of reviewing the various types of fractal 
mountains. The second type, the mountain range option, 
'Create a Mountain Range', creates a mountain range scene 
using the varicus types of fractal mountains. The first four 
options the Tutorial, Texture, Lighting Control ee 
Factor are for use with both a single mountain and a mountain 
range. The option Outline is for use only for creation ofa 


mountain range. Each of the options given above will be 


discussed in the order they appear. 


A. TUTORIAL 

The Tutorial option gives an online description of the 
different interactive options available for the creation of 
fractal mountains. The tutorial on the computer covers the 
same information as included in this chapter, but is less 
detail. The tutorial or online information is included so 
the user can use ISFM directly without reference to this 


chapter. 


B. TEXTURE 

Since the different types of mountains and mountain ranges 
vary in nature, fractal mountains need to be specified by a 
variety of textures and colors in order to achieve a realistic 
variety of fractal mountains. The option Texture gives five 
types of defined textures and associated colors. These five 
types are provided as options from a pop-up submenu of the 


Texture option. 
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The submenu is: 
Texture 
rugged granite 
rocky mountains 
alpine 
appalachian 
smooth rolling hills 
As mentioned in Chapter II.B.4, the ruggedness or 
smoothness of the texture is determined by a variable called 
texture. The variable texture determines the range permitted 
in the displacement of the midpoints. The larger the value 
of texture, the greater the displacement of the midpoints 
allowed and the more rugged the appearance of the fractal 
mountain. The smaller the value of texture, the smaller the 
displacement of the midpoints and the smoother the appearance 
of the fractal mountain. The orientation of the base 
triangles is also given for each option. These orientations 
are selected to match the type of mountain created. A further 
discussion of the orientation of the base triangles is found 
below under the create a mountain description. Each of the 
texture types, along with the color for the mountain, are 
discussed in the order given in the submenu and are all shown 
in the summary table, Table 6.1. Sample pictures are given 
in the results chapter (VII) for all the mountain types. 
1. Rugged Granite 
The rugged granite option yields a rock and boulder 


type texture, The texture variable is set to 15.5. In 


Chapter V, we state that the surface light properties consist 


61 


of ambient light, diffuse light and specular light. These 
three light types give the surface color for the mountain or 
mountain range. The typical values selected for the granite 


surface color are-ambient light red = 0.400, green = 0.415 and 


blue = 0.505. Diffuse light values are red = 0.405, green 
0.405 and blue = 0.500. Specular light values are red = 
0.464, green = 0.364 and blue = 0.564. Some specular 
component is included since a granite surface may be smooth 
and have directional light reflection. A foreground color is 
used to give the perception of depth. The color is given with 
red, green and blue components. The values of these three 
components range from O to 255. The color for granite is red 
= 175, green = 160 and blue = 178. 

The selection of the orientation of the base triangles 
was made in order to achieve a rugged outline for the 
mountains. In the single mountain case, the values of z depth 
range for all the mountains are from -900.0 (close to the 
eye) to -699.0 (far from the eye). The difference between the 
different textures appears in the y direction. For the rugged 
granite, the y values range from 200.0 to 550.0. For the 
rugged mountain range, the values of z depth range is -1020.0 
at y=0 to -700 at y=850. 

2. Rocky Mountains 

The rocky mountain option yields a rocky sharp 

texture. The texture variable is set to 12.5. The surface 


light properties consist of ambient light, diffuse light and 
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specular light. The typical values selected for the rocky 
mountain surface color are ambient light red = 0.565, green 
= 0.335 and blue = 0.035. Diffuse light values are red = 
0.585, green = 0.415 and blue = 0.325. For specular light all 
three components are set equal to 0.0. The foreground color 
for rocky mountain is red = 156, green = 102 and blue = 58. 

The selection of the orientation of the base triangles 
is make to achieve a rocky mountain. In the single mountain 
case, the values of the z depth range is again from -900.0 to 
-699.0. For the rocky mountain, the y values range from 200 
to 500. For the mountain range, the z values range from - 
1020.0 where y = 0.0 and go to -580.0 where y is 850.0. 

3. Alpine 

The alpine option yields a smoother mountain top 
texture than the rocky mountain option. The texture variable 
is set to 8.5. The surface light properties consist of 
ambient light, diffuse light and specular light. The typical 
values selected for the alpine surface color are ambient light 
red = 0.665, green = 0.285 and blue = 0.005. Diffuse light 
values are red = 0.725, green = 0.565 and blue = 0.235. 
Specular light is set equal to 0.0 for all three components. 
The foreground color for alpine is red = 201, green = 102 and 
blue = 58. 

The selection of the orientation of the base triangles 
is made to achieve an alpine mountain. In the single mountain 


case, the values of the z depth range from -900.0 to -699.0. 
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For the alpine mountain, the y values range from 200 to 450. 
For the mountain range, the z values range from -1020.0 where 
y= 0 and go to -250.0 where y is 850.0. 

4. Appalachian 

The appalachian option yields forested texture. The 
texture variable is set to 5.0. The surface light properties 
consist of ambient light, diffuse light and specular light. 
The typical values selected for the appalachian surface color 
are ambient light red = 0.015, green = 0.465 and blue = 0.025. 
Diffuse light values are red = 0.105, green = 0.545 and blue 
= 0.115. The specular light is set equal to 0.0 for all three 
components. The foreground color for appalachian is red = 61, 
green = 95 and blue = 36. 

The selection of the orientation of the base triangles 
is make to achieve an appalachian mountain. In the single 
mountain case, the values of the z depth range from -900.0 to 
-699.0. For the appalachian mountain, the y values range from 
200 to 400. For the mountain range, the z values range 
from -1020.0 where y = 0 and go to -60 where y is 850.0. 

5. Smooth Rolling Hills 

The smooth rolling hills option yields grass texture. 
The texture variable is set to 2.0. The surface light 
properties consist of ambient light, diffuse light and 
specular light. The typical values selected for the smooth 
rolling hills surface color are ambient light red = 0.025, 


green = 0.905 and blue = 0.035. Diffuse light values are red 
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— 0.245, green = 0.994 and blue = 0.155. The specular light 
is set equal to 0.0 for all three components. The foreground 
color for smooth rolling hills is red = 0, green = 95 and blue 
= 40. 

The selection of the orientation of the base triangles 
is make to achieve a rocky mountain. In the single mountain 
case the values of the z depth range from -900.0 to -699.0. 
For the smooth rolling hills the y values range from 200 to 
300. For the mountain range the z values range from -1020.0 
where y is 850.0. 


where y = O and go to -70 


Table 6.1 TEXTURE AND LIGHTING CONDITIONS 


Granite Rocky Alpine Appalachian Rolling Hills 

Texture 19.9 1225 8.5 5.0 2.0 
Ambient 

red 0.400 0.565 0.665 0.015 0.025 

green 0.415 02335 0.285 0.465 0.905 

blue 0.505 0.035 0.005 0.025 0.035 
Diffuse 

red 0.405 0.585 02725 0,105 0.245 

green 0.405 0.415 0.465 0.545 0.994 

blue 0.500 02325 0.235 0. 5 0.155 
Specular 

red 0.464 0.0 0.0 0.0 0.0 

green 0.364 0.0 0.0 0.0 0.0 

blue 0.564 0.0 0.0 0.0 0.0 
Foreground 

red 175 156 201 61 0 

green 160 102 102 95 95 

blue 178 58 58 36 40 
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C. LIGHTING CONTROL 
The lighting is composed of the surface light and color 
properties, the light source and the light level for the whole 
scene. In this section, the light source and the light model 
are discussed. The light source for the fractal neta is 
assumed to be the sun. As such, the light source is placed 
at infinity. Before creating the fractal mountain, pop-up 
menus give the option to set the time of day. 
The menu is specified as: 
Lighting Control 
Morning 
Noon | 
Evening 
Each of the above options are discussed separately and 
then a summary table of all the values are given in Table 6.2. 
1. Morning 
For the mountain or mountain range to be viewed in the 
morning light, the sun location must be in the east and at a 
low angle from the horizon. The vector x, y, z and w for the 
light source has values for morning light of x = -1.0, y = 
1.5, z= 0.5 and w = 0.0. With the variable w set equal to 
0.0, the light source is automatically set at infinity. The 
other three values determine the light direction. In the case 
of morning light the numbers indicate that the sun is coming 
in from the left side of the screen display at about a 45 


degree angle from the horizontal negative x axis. It is also 
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oriented about 40 degrees out of the xy plane. For 40 degrees 
out from the xy plane the lighting for the tilted base 
triangles is achieved. 

Two other variables determine the light source. These 
variables are the ambient light and the color of the light. 
The ambient light values for the morning scene is low because 
there is not much reflecting light in the morning. The values 


for ambient light are set to red = 0.2, green = 0.2 and blue 


0.1. The color for the light source is always the same no 
matter the time of day. The values are set to red = 1.0, 
green = 0.8 and blue = 0.0. 

The last component of the lighting is the light model. 
The properties in the lighting model are the scene ambient 
light and the viewer's location. Since the ambient light in 
the morning is low, the values are set to red = 0.2, green = 
0.2 and blue = 0.2. The viewer is placed at infinity, since 
this is the fastest lighting calculation mode. 

2. Noon 

The next option is the noon lighting option. The sun 
is placed directly behind the observer and at an angle of 60 
degrees from the horizontal. This direction of the sun is 
achieved by the values of x = 0.0, y = 2.0 and z = 1.0. Again 
w is set equal to 0.0, so the sun is located at infinity. 
The ambient light from the light source is considered to be 
higher at noon than in the morning. Values for ambient light 


are selected as red = 0.6, green = 0.5 and blue = 0.5. The 
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color of the light source is the same as for morning that is 
red = 1.0, green = 0.8 and blue = 0.0. 

The light of the scene for noon is the same in the 
morning except far the ambient light. The ambient light for 
the light source is larger at noon then in the morning. 
Therefore, the scene ambient light should also be larger. 
Values for ambient light of the scene are selected as red = 
0.6, green = 0.6 and blue = 0.6. 

3. Evening 

The last option is the evening light. The values of 
x, yand z are 1.5, 1.0 and 1.0, respectively. This results 
is about 30 degrees from the horizontal and about 30 degree 
angle out from the xy plane. The difference in the direction 
of the light from the morning compared to the evening is 
because of the tilt of the base triangles. Also in the 
morning the light comes from the negative x direction or east 
and in the evening the light comes from the positive x 
direction or west. 

The ambient light in the evening typically has a red 
glow, so the value of red is left at 0.6 as in the noon 
ambient light. The green and the blue values compared to the 
noon ambient light are decreased to a green value of 0.2 and 
a blue value of 0.1. 

Also, in the evening light was the amount of scene 
ambient light is changed. The evening ambient light model is 


more than the morning ambient light, but less than the noon 


68 


ambient light. The value for all three colors, red, green and 


blue, is 0.4. 


Table 6.2 LIGHT SOURCE AND SCENE LIGHT FOR LIGHTING CONTROL 
Light Properties Morning Noon Evening 


Light Source 


Position 
x SO 0.0 1.5 
y 1.5 2.0 120 
Z 0.5 1.0 13 
W 0.0 0.0 0.0 
Ambient Light 
red Orr 0.6 0.4 
green 0.2 0.4 0.2 
blue 0.1 0.4 0.1 
Color 
red 1:0 1:0 1.0 
green 0.8 0.8 0.8 
blue 0.0 0.0 0.0 
Scene Light 
Ambient Light 
red 0.2 0.5 03 
green 0.2 0.5 0.3 
blue 0.2 0.5 0.3 
Local Viewer 0.0 0.0 0.0 


The result of each of the three light sources and 
light models is that the triangles with different surface 
normals are each colored with a different shade. The 
different shades of color result help in achieving the depth 


perception of the fractal mountain. 
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DR STRUCTURAL FACTOR 
The structure factor option permits variation in the 
detail of the surface structure. The structure detail is 


determined by a submenu which gives the options: 


Structural Factor 
Medium 
Fine 
1. Medium 
When the structural factor is set to medium, the 
recursion stops when the length on an edge (using only the x 
and z coordinates of the two endpoints) is less than or equal 
to 16 pixels. If the length of the edge is larger than 16 
pixels for any of the three edges of the triangle, then 
another recursive Call is made for that triangle. If the 
length of all three edges is less than 16, then the triangle 
is colored, lighted and sent to the screen. 
2. Fine 
The other option for the structural factor is fine. 
The fine structure causes at least one more recursion level 
to be produced than the medium structural factor. The number 
of recursion levels is larger because instead of stopping at 
the length on an edge of 16 pixels as for the medium 
structure, the fine structure uses a stopping rule of an edge 


length of less than or equal to 9 pixels. The result of using 
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9 pixels for the edge length means a finer detail. The trade 
off is that it takes more time to generate the fractal 
mountain. 

The medium structure can be used for rugged to alpine 
mountains. The fine structure should be used for rolling 
hills. With the fine structure, the triangles used in the 
construction of the fractal mountain become less noticeable. 
The problem with using the fine structural factor is that it 
takes longer to generate the fractal mountain compared to the 


medium structure. 


E. CREATE A MOUNTAIN 

The create a mountain option is used for creating a single 
mountain. This option can also be used as a tutorial in the 
different types of fractal mountains available. In the 
creation of the mountain, only one base triangle is used. 
This base triangle is typically tilted out of the xz plane, 
so that the displacements in the y direction are seen more 
distinctly by the observer. The tilt is accomplished by 
placing the triangle on a plane going up in the y direction 
and going into the screen in the z direction. The amount of 
change needed in the y direction or in the z direction is 
determined by the type of mountain desired. 

The shape of the tilted base triangle as projected on the 
Xy plane is a factor in determining the type of fractal 


mountain created. With a base triangle having an acute apex 
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angle, such as triangle I in Figure 6.1, a sharp peaked 
mountain is created. If, on the other hand, the point B is 
changed in the y direction to B' with the same value of x and 
if the points A and C are kept constant, then the resulting 
triangle AB'C is closer to an equilateral ae 


triangle II). 





Figure 6.1 Types of Base Triangles 


An equilateral triangle gives a smaller or less sharp peak 
in comparison with the original triangle I. Further 
decreasing B' to B'' as shown in triangle III results in the 
triangle AB''C which is an obtuse triangle. The obtuse 
triangle results in mountains which have smooth or rounded 


apex angles or tops. 
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A combination of change in the y direction and changes in 
the z direction contributes to the type of texture desired in 
the fractal mountain and the realism of the fractal mountain. 
To expand the mountain to a more realistic view of nature, the 


mountain range option is included in ISFM. 


F. OUTLINE 

The outline section is used only with Create a mountain 
range. The outline allows the observer to create the contour 
of the mountain range. There is a submenu for drawing the 


outline which is: 


Mountain Outline 


create 
clear and reset 


The main function of the outline section is to place the 
mountain range wherever the observer wants. The create 
submenu gives that capability. 

1. Create 

The create option gives a clear screen with the 
initial point in the center of the window. To create the 
outline of the mountain range move the cursor along and press 
the left mouse button. If the cursor moves less than 5 pixels 
in either the x or y direction than the red line will not be 


drawn from the previous point to the cursor. The value of 5 
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is a noise value so that not every slight movement of the 
cursor is reported. Once the left mouse button is pressed the 
beginning and ending values of the red line is stored and is 
then shown as a-black line on the screen. The red line now 
begins at a new location, i.e., the last x and y values. The 
process of moving the cursor and pressing the left mouse 
button can continue across the window until the outline is 
made to fill the whole width of the window. The outline must 
be a continuous line starting from the left and running to the 
right or starting from the right and running to the left. 
Since the location of the cursor begins in the center 
of the window and the outline must be continuous, the 
capability of replacing the original starting point must be 
available. To replace the starting point or to create a new 
location for the x and y values the middle mouse button is 
used. To achieve a new location with the middle mouse button 
the cursor is moved. This stretches the red line, to the 
point where the outline is to start. Once the cursor is 
placed in the desired location, the middle mouse button is 
depresses and the new location is set. The result is that the 
red line now starts at this new location. From this new 
location the rest of the mountain outline can be drawn. 
In order to be able to change the outline after having 
created all or part of the outline, the right mouse button is 
used. That is, every time the right mouse button is depressed 


the last point drawn is erased. This erasure of the last 
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point can be continued until there are all points are erased. 
This feature permits changing the outline after it has been 
created. 

After the outline has the desired shape for the 
mountain range, there is the option of exiting the outline. 
To exit press 'e' on the keyboard. Once the 'e' is pressed 
the given outline will stay up on the screen as a black 
contour. 

If one wants to change the outline, the create option 
in the submenu Mountain Outline can be chosen again. If the 
create option is chosen again the cursor will start at the 
last point from the previous outline. If on the other hand, 
one wants to create a completely new outline press the other 
option clear and reset. 

2. Clear and Reset 

The clear and reset option is used to clear the whole 
window. In this way one obtains an empty window in which a 
new outline can be created. This option also resets the 
values in the arrays to zero to start with a reinitialized 
array. This allows one to create a completely new outline for 


the mountain range. 


G. CREATE A MOUNTAIN RANGE 
The create a mountain range option is an extension of the 
create a mountain option. Instead of one base triangle, the 


mountain range is created with multiple base triangles. The 
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multiple base triangles each have a change in the y direction 
and a change in the z direction, similar to what was done in 
the case of one base triangle. The combination of the two 
changes contributes to the whole surface or texture of the 
mountain range. 

The area in which the mountain range is to be placed is 
covered with triangles. That is, the triangles make a pattern 
on the screen display or as a projection on the xy plane, as 
in Figure 6.2. The triangles are numbered in the order that 


they create the fractal scene. 





Figure 6.2 Layout of Base Triangles for Mountain Range 


In Figure 6.2, the process of creating the mountain range 
begins in the bottom left hand corner and moves across the row 


alternating the triangles from one vertex pointing up to one 
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vertex pointing down. The values assigned to each of the 
endpoints of adjacent triangles are the same in order to avoid 
the appearance of the gapping problen. Each vertex of the 
base triangle pattern is a function of x, y and z coordinates 
as defined in our usual coordinate systen. Each of the 
triangles 1 through 5, and 7 through 11 are all of the same 
size. The y coordinates can be varied so that one has 
triangles with apex angles from obtuse to almost equilateral 
angles. The size of the triangles is selected to match the 
type of texture desired. Triangles 6 and 12 are placed on the 
sides and configured in order to have the size of all the 
triangles be about the same. The size is important because 
the number of recursion levels in each triangle depends on the 
size of the original base triangle. If the number of 
recursion levels in each base triangle is not approximately 
the same, then another form of the gapping problem occurs. 
This gapping problem occurs when considering for the adjacent 
triangles the midpoint for one triangle is displaced but not 
for the other adjacent triangle if the number of recursion 
levels are different. In order to assure that each base 
triangle uses the same number of recursion levels, the number 
of recursion levels for the first base triangle is determined 
based on the size of the smallest triangle desired. For all 
succeeding base triangles, the stopping criterion is based on 
the same number of recursion levels instead of the size of the 


smallest triangles. It is for this reason that the pattern 
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of covering base triangles has been composed of essentially 
equal size triangles. 

The size of the screen display and the magnitude of the 
angle of the obtuse triangle used determines the number of 
total base triangles used to create the entire mountain scene. 
The pattern of twelve triangles as shown in Figure 6.2 is 
repeated vertically until the entire screen display is 
covered. 

One of the important elements in creating a mountain range 
is the depth perception. Depth perception is achieved by 
using different z values, as described in Chapter V. To 
provide depth perception to the fractal mountain, the values 
of z depend on the values of x and y. As the y values 
increase, the z values also increase. If y was the only 
coordinate upon which z depended, then depth perception would 
be directly into the screen. Instead, the z values also 
depend on the x coordinate. As the x values increase, the z 
coordinates also increase. The amount of change in the z 
coordinate is affected more by the y coordinate than by the 
x coordinate. This enhances the natural look of the scene. 
With the z coordinate dependent on both x and y coordinates, 
one has a sense of depth perception and of perspective which 
is into the screen and towards the right. 

The dependence of the z coordinates on both the x and y 


coordinates helps create a realistic orientation for the 


78 





surface of the fractal mountain or mountain range. Variations 
in the orientation yield different viewing directions. 

As explained in the next section color variations and flat 
surfaces are also introduced in the mountain scene in order 


to create a further perception of depth. 
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VII. RESULTS SHOWING FRACTAL MOUNTAINS 


ISFM gives E capability of creating various types of 
mountains. In order to give depth to the pictures, a 
variation from light to dark color blue was put into the sky. 
Also a solid foreground polygon is created to give the depth 
in the pictures. There are a couple other aspects that help 
the realistic look of the fractal mountains. These aspects 
are weighted displacement of the midpoints and sophisticated 
lighting. The weighted displacement of the midpoints allows 
for less displacement of the midpoint as the triangle size 
gets smaller. The sophisticated lighting allows the mountains 
to be see in the morning, at noon and in the evening. For 
creating the mountains with the options of texture, lighting 
control and structural factor as well as the aspects just 
mentioned, one can create a single mountain or a mountain 


range. 


The results of ISFM of each of the sections in Chapter VI 
is shown by pictures in the following sections. The first 
option in ISFM is the texture section. The different types 
of textures are used to show some of the other options 
available in ISFM. To show the textures from roughness to 


smoothness, the single mountain case is used. In the single 
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mountain case, there is a solid foreground corresponding to 
ene color of the type of mountain. The sky color is shaded 


from a light blue to a shade of dark blue with increasing 


> 





Figure 7.1 Rugged Granite Mountain 


height to give the observer the perception of depth. Figure 
7.1 gives an example of the most rugged type of texture, that 
of rugged granite mountain. 

The next option of texture, the next step down in 
roughness, is the texture called rocky mountain. The example 
of a rocky mountain is shown in Figure 7.2. In comparison to 
the granite mountain, the texture is different and the 


mountain is not as high. 
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Figure 7.2 Rocky Mountain 


To show the lighting control options, the texture of 
alpine mountain has been used as an example. The lighting 
options are morning, noon and evening. The morning light on 
an alpine mountain is shown in Figure 7.3. As can be seen the 
light is coming in from the left side of the observer. 
Figures 7.4 and 7.5 are examples of the light of noon and 
evening, respectively, on an alpine mountain. Notice the 
light on the mountain at noon comes from behind the observer 
and the light at the evening comes from the right side of the 


observer. 
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Figure 7.3 Alpine Mountain - Morning Light 





Figure 7.4 Alpine Mountain - Noon Lighting 
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Figure 7.5 Alpine Mountain - Evening Light 


The next option is the structural factor. In order 
give an example of the two structural factors, the texture 
the appalachian mountain is used. The structural factor 
medium is shown in Figure 7.6 and the structural factor 


fine is shown in Figure 7.7. 


to 


of 


of 


of 


The structure of Figure 7.7 is finer than the structure 


of Figure 7.6. The finer structure is due to the existence 


of at least one more recursion level in the mountain formation 


than with the medium structure option. 
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Figure 7.6 Appalachian Mountain - Medium 





Figure 7.7 Appalachian Mountain - Fine 
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An option available only in create a mountain is the 
magnification factor. This magnification factor allows one 
to see the structure of the mountain in more detail. 

The last texture type of smooth asii hills is used to 
demonstrate how the magnification looks, Figures 7.8, 7.09, 
7.10 and 7.11 have the magnification levels of 1.0, 2.0, 3.0 
and 4.0, respectively. These magnifications give the observer 
a good idea of what the structure actually looks like down to 


the pixel size. 





Figure 7.8 Rolling Hills Mag = 1 
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Figure 7.9 Rolling Hills Mag = 2 





Figure 7.10 Rolling Hills Mag = 3 


87 





— _ A A A AA TRA nn DEE Rn ut aac 





Figure 7.11 Rolling Hills Mag = 4 


The above figures give all the options available for a 
single mountain. The other two options that of outline and 
create a mountain range are shown in a series of pictures of 
the mountain range given for different outlines. The outline 
of each mountain range is indicative of the type of texture 
of the mountain range. The texture again is given from the 
most rough to the smoothest. Each of the options of lighting 
control and structural factor is incorporated in one or more 
of the examples. Figure 7.12 is an example of a rugged 
granite mountain range with a morning light and a structural 


factor of medium. 
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Figure 7.12 Rugged Granite Mountain Range - Morning Light 
- Medium 


The next example, Figure 7.13, is of a rocky mountain range 
with a noon lighting an a medium structural factor. 

Figure 7.14 is an example of an alpine mountain range with 
evening light and a medium structure factor. 

To show an example of an appalachian mountain, the options 
of morning light and fine structure factor are shown in Figure 
2275: 

The last mountain range texture is the smooth rolling 
hills, as seen in Figure 7.16. The smooth rolling hills 
picture has the options of an evening light and a fine 


structural factor. 
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Figure 7.13 Rocky Mountain Range - Noon Lighting - Medium 





Figure 7.14 Alpine Range - Evening Light - Medium 
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Figure 7.15 Appalachian Range - Morning Light - Fine 





Figure 7.16 Rolling Hills - Evening Light - Fine 


ot: 





The above pictures Figures 7.1- 7.16 give the available 
options for both creating a mountain and creating a mountain 
range. These figures are a representation of the different 
possibilities available while using ISFM. The results display 
a good deal of realism and natural effects typical of mountain 


scenery. 
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VIII LIMITATIONS, FUTURE DIRECTIONS AND CONCLUSION 


The Interactive System for Fractal Mountain (ISFM) is 
user friendly and facilitates the creation of Fractal 
Mountains. It has the capability of creating various type of 
mountains. In order to give depth to the pictures, a 
variation from light to dark color blue was imparted to the 
sky. Also a solid foreground polygon is created to lend the 
depth to the pictures. The use of sophisticated lighting 
routines and depth perception is enhanced by coloring and flat 
Space segments in the mountains and the mountain range. 

The Interactive System for Fractal Mountain (ISFM) can be 
expanded to include seasonal effects for the mountain ranges. 
The seasonal effects can be accomplished by varying the basic 
grid pattern for the mountain range. Other things that could 
be done with the system, are to introduce other natural 
aspects such as lakes and rivers within the scene. Along with 
including lakes and rivers, a whole mountain scene could be 
included by a progression of different textures across the 
basic grid pattern. Different textures within a mountain 
scene would allow the creation of rolling hills and rocky 
mountains to be in the same picture. To permit the different 
textures within the same picture to be chosen by the user, the 


user would need to be able to work on different grid sections 
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of the mountains independently. Working on different sections 
could also allow the possibility of color variation across the 
different sections of the picture. This would allow, for 
example, the creation of snow-capped mountains by the user. 

Two other possible extensions to the interactive system 
concern the outline and the foreground. The outline should 
match the type of mountain being created. This could be done 
by randomly displacing the points of the given outline in the 
same manner as used in the displacement for the mountain 
texture. Another useful addition to the interactive system 
would be the ability to vary the size and shape of the 
foreground contour. This variation would give different depth 
perceptions to the user. 

Not only can this interactive system be expanded to 
include more capabilities for the fractal mountains, but also 
the system could be used to create other natural objects like 
clouds, trees and ocean surfaces. These other natural objects 
could use some of the same techniques used in creating fractal 
mountains. One of the possible applications could be based 
on different base shapes. In this study, a triangle is used 
as the base shape. To create other natural objects, shapes 
such as squares, rectangles or any other polygon could be 
used. 

The Interactive System for Fractal Mountain (ISFM) is a 
straightforward method of developing fractal mountains. ISFM 


gives the capability of creating various types of mountains. 
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ISFM is only the first part of a project whose ultimate goal 
is the complete description of the techniques necessary for 
utilizing fractals for the generation of realistic texturing 


for computer generated images. 
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APPENDIX : C PROGRAMS 


PROGRAM : MNTMAIN.C 


/* This is program mntmain.c 
sets up the pop-up menus for Fractal mountains */ 


#include "gl.h" 
#include "device.h" 
#include "list.h" 
#include "math.h" 
finclude "stdio.h" 


define EXIT 99 


/* integer variables for menus - main, mountain outline, 
mountain texture, lighting control, structural factor, 
magnification and exit ar 


int mainmenu; 
int mntmenu; 
int txtrmenu; 
int lsmenu; 
int strtmenu; 
int exmenu; 
int crtmenu; 


/* integer variables denoting the current texture, structural 
factor and light position */ 


int ceurlsor: 
int eurtxtr: 
aäne=strtfer; 


/* variables denoting the current magnification factor 
texture and single mountain or mountain range S 
int imag; 
float texture; 
int m; 


/* Counters */ 
ine IC Tris 
inesi. 


/* pointer variables for link list which stores the vertices */ 
LINK first, frst,tail; 
float SR; 


/* array to store the vertex values 2 
float ASIS 

/* array for the random numbers */ 
double RAND[500]; 
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/* arrays for storing the outline and integers for total amount 
of points in each array ne 
Pleat prts[ 10007713]: 
anemotrcps;: 
eisat alpts [1000] [3]; 
Intsntot; 


tr 


long op[1000]; rs Thereperatien to perform for the outline */ 


/* flags for mountain outline, number of levels, logical variable 
if 1 or more base triangles have been created, default outline used 
and how many triangle patterns used in creating mountain range */ 
ine flg; 
int nlev,ct; 
nt Ot © Lg; 
float squfac; 


/ EOmcedon to Set the Tight source at morning, noon or evening */ 
1sp (n) 
( 
Gurlsor =r; 
return 3; 


} 


/* £unction to set the texture type rugged granite to smooth rolling hills */ 
txtrmnt (n) 


ine n; 

{ 

eurextr = N; 
geturn 2: 


) 


/* function to set Structural factor medium or fine */ 
stertmnt (n) 


ine n: 

{ 

strtfct = n; 
return 4; 


) 


unction to set the magnification of the single mountain */ 
crtmnt (n) 


int n: 

{ 

imag = n; 
return 5; 


) 
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main () 

{ 

/* value for reading from the queue */ 
short value; 

inei: 


/* integer value for pop-up menus a 
static int pupval = 1; 


/* create and open Fractal Mountains Window */ 
prefsize(1050,850); 


winopen ("Fractal Mountains"); 
wintitle("Fractal Mountains"); 
winconstraints(); 


/* Place in the modes singlebuffer and RGB 
and set the window depth and screen coordinates */ 


singlebuffer (); 

RGBmode (); 

viewport (0,1050,0,850); 
lsetdepth(0x0000,0x7fffff); 

ortho (0.0; 1050.0,0.07850.0, 102370 070 
geonfiaq(), 

qdevice (REDRAW) ; 

qdevice (MENUBUTTON) ; 


/* Build the random number table */ 
rand table gen(); 


RGBcolor (0,175,255); /* color blue*/ 
clear (); 
/* create the intro window with text “fl 


RGBcolor (0,070) 

cmov2i({275,600):; 

charstr ("Welcome to Interactive System for Fractal Mountains(ISFM)"); 
cmov2i (250,450); 

charstr("If this is your first time please see the Tutorial section"); 
emov21.(250,400): 

charstr ("for a description of the options and defaults"); 
cmov22.(2007150); 

charstr("Program designed and written by Patricia Collins-Hunt"); 


/* Create the Second level exit menu */ 
exmenu = defpup ("Exit %t |Yes I am sure $x99|No don't do it"); 


/* Create the menu for Position of the Light Source Bf 


lsmenu = defpup("Time of Day %t %F |morning %x3l|noon %x32 ",1sp); 
addtopup(lsmenu, "evening %x33") ; 
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/* Create the menu for mountain outline = 
mntmenu = defpup ("Mountain Outline %t"); 

addtopup (mntmenu, "create %x61"); 

addtopup (mntmenu, "clear screen %x62"); 


/* Create the menu for Magnification = 


ertmenus=zdefpupi(7Magnifsceastien st %E|IT.O %x1 |2.0 %x2 [3.0 %x3",cremnt): 
addtopup (crtmenu, "4.0 %x4"); 


/* Create the menu for mountain texture ay 


txtrmenu = defpup ("Mountain Texture %t %F | rugged-granite %x21",txtrmnt); 
addtopup (txtrmenu, "rocky mountains %x22 |alpine %x23 |appalachian %x24" ); 
addtopup(txtrmenu, "smecth rolling hills %x25" ); 


/* Create the menu for Structural Factor “7 
strtmenu = defpup ("Structural Factor %t %F | medium %x16",strtmnt); 
addtopup (strtmenu, "fine %x9" ); 


/* Create the main menu for Fractal Mountains = 
mainmenu = defpup ("Fractal Mountains %t |Tutorial %x1"); 
addtopup (mainmenu, "Texture %m",txtrmenu) ; 
addtopup (mainmenu, "Lighting Control %m",lsmenu) ; 
addtopup (mainmenu, "Structural Factor %m", strtmenu) ; 
addtopup (mainmenu, "Create a Mountain %m ", crtmenu); 
addtopup (mainmenu, "Outline %m ",mntmenu) ; 
addtopup (mainmenu, "Create a Mountain Range %x7"); 
addtopup (mainmenu, "clear and reset %x8jexit %x97"); 


/* initialize the variables denoting the current options */ 
curtxtr = Q0; 
cuüurlisor = 0; 
Strt fct 0; 


/* Go into loop for pop-up menus until exit is chosen */ 
flg = 0; 
while (pupval != EXIT) 
( 
switch(qread (&value) ) 
( 
case REDRAW: 
reshapeviewport (); 
break; 


case MENUBUTTON: 
if (value == 1) 
( 

pupval = dopup (mainmenu); 

switch(pupval) 

( 

case 1: 7A for Tutorial */ 
etutorl(); 
break; 


99 


case o: /* create a mountain nd 
m=1; 

sampmnt (); 

break; 
case 7: /* create a mountain range */ 
m=2; 

sampmnt (); 

break; 
case 8: /* clear and reset option */ 
RGBE919r210,175,255), 

clear (); 

break; 


case 61: /* for creating outline for mountain range */ 
drawmnt () ; 

break; 

case 62: /* for clear and reset outline option a 
RGBcelor (0, 17525255) 

clear (); 


for (i=0;i<=1000;1++) 
( 
alpts[(i][0]=0.0; 
alpts[i])](1]=0.0; 
alpts[i] [2]=0.0; 


pts [1 10) =0707 
pts[i]([(1]=0.0; 
pts[i] [2]=070; 


flg = 0; 
) 
break; 
case 97: /* for exit option = 
pupval = dopup (exmenu); 
break; 


} 

} 
break; 
default: 
break; 


) 
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PROGRAM 


jr 


MOUNTDYN.C 


nis program initializes the system to light and color the fractal 
mountain or the mountain range */ 


/* The variables are the following 
texture-gives the value of the roughness of the mountain 


MEA 


counter for how many times went through gendyn 


Tri - how many triangles are created 


first 
st 
AS[10 


, frst,tail - pointers to addresses of the array of 
oring the vertices 
]- stores the three endpoints of the triangles and puts 


them into dynamic storage 


ysno- 


variable to send triangle to screen 


cüreti value Of Current texture of options 


curls 


or - value of current light position of options 


m - determines if single mountain or mountain range 


squfa 
mo 
PITE? 


#include 
#include 
#include 
#include 
#include 
#include 


#define x 
#define y 
#define z 


/* Global 


extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 


c - how many sets of 12 triangles in pattern for 
untain range 
‚P3 - endpoints of the triangles */ 


<gi- -h> 

<device.h> 

<math.h> /* Standard UNIX include file for math library */ 
A 

"staro.n"” 

ita hede£s .h" 


0 
1 
2 


Tables */ 


float texture; 
meee, Tri: 

DENK first, frst,tail:; 
loarz25[10]; 

float ysno; 
intsceurtxtr; 
intseurlsor:; 

int m, 

float squfac; 

float magfac; 


/* BEGIN SAMPMNT */ 


sampmnt () 
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/* Local Variables. 
float P131, P2 IL 2317312 
Colorindex wmask; 


/*Enable alí bitplanes for writing*/ 
wmask=(1<<getplanes ()) -1; 
writemask (wmask) ; 


wintitle ("Mountain Scene”); 
mmode (MVIEWING) ; 

RGBmode () ; 

singlebufter (); 

drawmode (NORMALDRAN) ; 
gconfig(), 


/* Fractalize until desired precision */ 
IC = 0; 
Tri = 0; 


/* initialize the lighting values */ 
initialmat (); 
initialls(); 
inicral im); 


zbuffer (TRUE); 
zclear(); 


/* determine the type of texture */ 
if (eurtxtr =E0) curexerie 23: 


switch (curtxtr) 
{ 
case 21: /* rugged granite */ 
Imbind (MATERIAL, MATGRANITE) ; 
texture = 15.5; 
squfac = 2.0; 
break; 


case 22: /* rocky mountains */ 
1mbind (MATERIAL,MATROCKS) ; 
texture = 12.5; 
squfac = 2.0; 
break; 


case 23: /* alpine mountains */ 
lmbind (MATERIAL, MATSUNSET) ; 
texture = 8.5; 
squfac = 3.0; 
break; 
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case 24: Ikt appalactan mountains) */ 
lmbind (MATERIAL, MATFOREST) ; 
texture = 5.0; 
squfac = 4.0; 
break; 


case 25: /* smooth rolling hills  */ 
lmbind (MATERIAL, MATGRASS) ; ~ 
texture = 2.0; 
squfac = 5.0; 
break; 


) 
/* determine location of the light source */ 


if (curlsor == 0) curlsor = 31; 
switch(curlsor) 
( 
case 31: /* morning sun */ 
lImbind (LMODEL, LMODM) ; 
lmbind (LIGHTO, LSOURCM) ; 
break; 


case 32: /* noon sun */ 
imbind (LMODEL, LMODN) ; 
lmbind (LIGHTO, LSOURCN) ; 
break; 


case 33: /* evening sun */ 
imbind (LMODEL, LMODE) ; 
Lmbind (LIGHT2, LSOURCE) ; 
break; 


) 


/* detrmine if a mountain or a mountain range is desired */ 
Switch (m) 


{ 

case 1: /* for creating a mountain Di 
mountcal(); 
break; 


case 2: /* for creating a mountain range */ 
interpmnt (); 


basetri(); 
break; 


zbuffer (FALSE); 
wintitle("Fractal Mountains"); 


/* END mountdyn */ 
} 
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PROGRAM 


MOUNTCAL.C 


pe This program initializes the basetriangle for creating a mountain */ 


/* The external variables are the following: 


EOS 


- counters to determine number of base triangles 


Tri - number of triangles created 

nlev - total number of recursive levels 

ct - counter for determining more than one base triangle 
curtxtr - value from options of current texture 

imag - amount of magnification 


a7 


finclude <math.h> /* Standard UNIX include file for math library */ 
#include "list.h" 
finclude "stdio.h" 


define x O 
define y 1 
define z 2 


/* Global 


extern 
extern 
extern 
extern 
extern 


Tables */ 


ine TE Te, Er. 

LINK first;frst; tail; 
int nlev ct; 
intreurexrter; 

int imag; 


mountcal () 


{ 


/* Local Variables */ 


/* 


IR 


/* 


variables for the endpoints of the triangles 27 
float Pins] EZO I PITE 

int ilev; 

variables for foreground polygon A 


float £rdl[3],£rdr [377 Er) 37 Fra coca E 
float forg2[37,£09r93 [37], 20293213], forged | 32] orgs] forge eE 


variables for sky polygon us 
float skal[3}, skdze[ 3), skul 317, skuz | 37 
flost=eup[3] can (31 


initialize values for foreground polygon */ 
£frdl (000; 

frdl(1]=0.0; 

£frdl [Zee -30070: 


frdar[0]=1050.0:; 


frdr [1] 00: 
£rdr [222 300.8: 
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EXUL(0]+0-0:; 
Pu] 23000; 
pruiz]="=300.0> 


fee | O}—1050.0; 
feur (1 p=355.0; 
frurf2)]=°>-300.0; í 


forgi [0]=85.0; 
forgl[1]=265.0; 
forgl[2]= =300. 0; 


£9292[0)=200.0; 
forg2[1y=257.0; 
forg2[2J=: -300.0; 


Gongs [(Ojeso5 .0; 
forg3[1)=240.0; 
forg3[2]= -300.0; 


forg3a[(0)=470.0; 
forg3a(1)=297.0; 
forgsal2)= 30070; 
forg4[0]=600.0; 
forg4[1)]=320.0; 
forg4 [2 j= -300.0; 


forg5 [0 3745.0; 
forg5(1)=340.0; 
foegows) = —300.0; 


£org6 ( 0y=915.0; 
forg6 [ 11=360.0; 
forgot |= -300.0; 


J Create @@he initiating structure for the mountain */ 
switch (curtxtr) 
( 
case 211 /* rugged granite case */ 
Prix tet 30020: 
Pl{y) = 200.0; 
Parc) = -900:0; 


Peg) 450.0; 
Ez]? 55050; 
P2[z] = -699.0; 


P3[x] = 600.0; 
P3[y] = 300.0; 
P3[z] = -800.0; 


REBECLOT (136, 143,128); 
break; 
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case 22: 
Pl{x] 
Ey 
P1{z] 


P2 [x] 


PZ [y] 
F2[z] 


P3[x] 
P3[y] 
Balz] 


/* rocky mountain case = 7 
30070; 
20.020: 


-900.0; 


450.0; 
500 07 
-699.0; 


60010; 
300 FUE 
=80 0 707 


RGBcolor(75,42,19)7 


break; 


case 23: 
P1{[x] 
ELY] 
pliz] 


P2[x) 
P2 [y] 
P2[z] 


P3[x] 


ES 
Paiz] 


/* alpine case */ 
300507 
200.0; 


=9007 0, 


450707 
45020; 
-699.0; 


600.0: 
300.0; 
-800.0; 


RGBCOLOr (159 8870F 


break; 


case 24: 
Pix] 
Pity] 
Pl[z] 


P2[x] 
P2[y] 
P2[z] 


P3[x] 


ES [y] 
PS{z] 


/* appalachian case */ 
30030; 
200207; 


-900.0; 


4500; 
400.0; 
-699.0; 


600.0; 
300. 0; 
-800.0; 


RGBcolor (61, 1195-36); 


break; 


case 25: 
P1{x] 
Ey) 
PLIZ] 


/* smooth rolling hills case 
3007 07 
20007 


-900.0; 
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zu 


P2 x} 
EA 
E Zale 


PS i) 
Pa [y] 
Esz] 


RGBcolor (0, 95,40); 


break; 


ilev = 
Ti. => 
Ic = 0 


CE 0; 


O; 


Tri = 0; 


45020: 
300.0: 


693.0; 


600.0; 
250.0: 
-800.0; 


/* initialize counters and flag i 


e rea cemthe foreground structure = 


bgnpolygon(); 
vor tt ral) 
Vat Cl ru) 7 
Mor (forgi) 
yoe Eorg2); 
v3el(torg>3r. 
v2:(20:933); 
v3f(forg4); 
wstiforg)); 
varı(forg6); 
YItıerur); 
VSE (Era) > 

endpolygon(); 


J simnveraalize the sky structure */ 


skdl[0] 
skdl[1] 
skdl[2] 


skdr [0] 
skdr[1] 
skar| 2] 


skul[0] 
skul[l] 
skul[2} 


skur[0] 
skur[1] 
skur [2] 


020: 
20020; 
— 00: 


POSO: 
20070: 
-100.0; 


DIO: 
850.0: 
-100.0; 


1050.0; 
295097 07 
z100. 0; 
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cup [0] =2020 
cup[1] = 
cup[2] = 


e O 

O U 

Se GO “e 
a 


edn | 0]7 =20930: 
cani] 0.88; 
can] etico; 


/* create the sky structure */ 


bgnpolygon () ; 
c3f (cdn); 
v3f (skdl); 
v3f (skdr); 
C3 (cup): 
v3f (skur); 
v3f(skul); 

endpolygon(); 


/* create the fractal mountain * / 
mountain generate (P1,P2,P3,ilev) ; 
listlite (first); 

system("mag imag"); 
ct =i l; 
return; 


/* END mountcal */ 
) 
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PROGRAM : RANDTBL.C 


ER Tìis program tests the random number generation for the 
mountain fractal gaussian distribution */ 


#include <math.h> /* Standard UNIX include file for math library */ 


/* external global variables */ 


extern double RAND[500]; 


/* BEGIN RAND_TABLE_GEN */ 
rand tabrie gen) 


{ 


/* Local Variables */ 


ine Lerch: 

double UNF1, UNF2; 
double range,pi; 
Inet Eactor; 


peen. 1415926535; 
/* Determine the range for the random numbers of UNIX UCB */ 
range = 2; 
for (J=1; J<=14; J++) 
i range = range * 2; 
ee =i range = 1; 
/* Set the random number generator seed */ 
srand (475836); 


/* Create a Table for 500 entries */ 


form T=0; I<. 5007 Im I + 2) 
{ 


/* Get a uniform random number through the Unix C subroutine */ 


UNEL 
UNF 2 


N 


Sana); 
rand (); 
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/* Normalize the uniform random number to the intervalzror eee 


UNEL 
UNF 2 


UNF1 / range; 
UNF2 / range; 


/* Mold the uniform random variable into the approximate normal 
distribution */ 


factor = 1.0; 

if (log(UNF1) < 0.0) factor = =R; 

RAND[I] = sqrt (factor * (2.0 * log(UNF1))) * 
cos ((2*pi*UNF2)); 

RAND[I] = RAND[I) * factor; 

factor = 1.0; 

if (log(UNE2) < 0,0), factor = -1.0: 

RAND [I+1] = sqrt (factor * (2.0 * Irog (UNE TIDE: 
sin ((2*pi*UNF2)); 

RAND [(1+1] = RAND(1I+1] * factor; 

) 


return; 


1:0 


PROGRAM. 3 LIST.H 


/* This file is used by other programs in creating the dynamic 
variable. In this program we create the linked list to be store 
the points*/ 


ù 


#define NULL 0 
typedef float AP; 


struct jinkeder:st 
{ 
float AP; 
struct iinked list *prev; 
Seructlinkeg list *next; 
ir 


typedef struct linked list ELEMENT; 
typedef ELEMENT *LINK; 
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PROGRAM; “Listes 


/* This routine will pick up the stored values and then 
light the polygons */ 


/* The variables are the following: 
IC, Il are counters for number of times gendyn and listlite used 
frst,tail,lst - pointers for addresses in dynamic storage 
AS[10] - array for storing endpoints of the triangles in 
dynamic storage 
ay 
#include "list.h" 
#include "stdio.h" 
#include <math.h> 


extern int IC. IE 
extern LINK frst,tail; 
extern float AS[10]; 


listlite (head) 
LINK head; 


{ 
LINK lst; 
a ales aie 


IC=IC+1; 

lst = head; 

if (IC==1) frst = head; 
if (IC==1) lst = tail; 


/* retrieve the three coordinates of each point anda flag = 
while (lst != NULL) 
{ 

for (i=9;i>=0;i=i-1) 
( 

AS[i] = lst->AP; 

free (lst); 

lst = lst-> prev; 
) 


/* send endpoints to litemnt to light the triangles and send 
to the screen e 
litemnt (); 
) 
zeturn, 


} 
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PROGRAM 


BASETRI:C 


/* This routine will pick up the stored values and then 


light the polygons */ 


#define x O 
#define y 1 
#define z 2 


finclude "gl.hn" 
#include "device.h" 
#include "list.h" 
#include "stdio.h" 
#include "lightdefs.h" 
#include "math.h" 


/* the variables are the following: 
pts and alpts - arrays for values of a given outline 
totopgs - total # of points in alpts array 
ntot - total # of points in plts array 
RAND - array for the random variabales 
texture - value for each of the different textures 
first, tail and frst - pointers to addresses of the array for 
storing the vertices 
Il,IC - counters for new base triangles 
Tri - counts # of triangles created 
nlev - total # of recursion levels 
ct - flag for determining if one or more base triangles created 
squfac - how many time triangle pattern is repeated for 
mountain range 
curtxtr - current texture value from menu 
otflg - for default outline 
ny 
extern float pts[1000] [3]; 
extern int totops; 
extern double RAND[500]; 
extern float texture; 
extern LINK first,tail, frst; 
esterni int IL, IC,Tri:; 
extern float alpts[1000] [3]; 
extern int ntot: 
extern int nlev,ct; 
extern float squfac; 
extern int Cuüurtxtr:? 
extern int otflg; 
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Basetrı() 


| 


/* the variables are the following 


/* 


ilev - the current recursion level 

cont - vector of points to create each vertex of a polygon 
1, J, k - loop counters for pattern of mountain range 

fg - counter for the current base triangle in the pattern 
maxtri - maximum number of base triangles created 


vupr, vupl, vdownr, vdownl, vnxtr, vnxtl - arrays for creating 
the sky polygon 
cvd, cvup ~- color arrays for the sky polygon 


ts - variable to make sure outline starts at 0 and ends at 1023 

vn - 1 or -1 for creating each series of 12 triangles 

yfac - a factor used in the pattern for the mountain range 

P1, P2, P3 - values of endpoints of the triangle 

TP - temporary array holding values of the endpoints 

loc - variable to determine location in array RAND 

frmdl, frmdr, frmul, frmur = arrays for greatingrforegroundpelvaer 

frml - frm6 - arrays for creating foreground polygon 

znn, zp, zn, zpp - variables for changing z values for each texture 
=) 


int ilev; 

float cont[3]: 

int i; k: 

int fg mart ros 

float vupr [3]; vupl (3 leedornc (slp vacunl ier 
float. vnxexr lor vonke MoE 

float teva la], cruel 

float ts, vn; 

float yfac; 

float P1(3),P27 9, P3 [2], ER 3]: 

int ae: 

float 7frmall[l], frmarl 3), Zrmul 2 ana 
float £rml1 [3], £rm2([3], £xrm3 [3], txm4 (3), tres see ee 
float zan, Zp; zn, Zp; 


winconstraints(); 
initialize values for the sky polygon a 
vupr[0]=1050.0; 
vupr[1]=1050.0; 
vupr[2]= -1023.0; 
vupl [0] =020; 


vupll[1]=205070; 
vupl[2]= -1023.0; 
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/* 


vdownr[0]=1050.0; 
vdownr[1)]=0.0; 
zdennr[2]= =1023.0* 


vdownl[0]=0.0; 
vdownl(1]=0.0:; 
vdownl(2]= -1023.0; 


vnxtr([0])]=1050.0; 
vnxtr[1]=10.0; 
vnxtr [2)]= -1023.0; 


vnxt1(0]=0.0; 
vnel [1]=100.0; 
wheelie 2)— -=1023.0; 


cvd[0] 2.0; 
cual] =40.88; 
erd(2]j 201.0; 
cvup{0) = 0.0; 
erupii]) = 0.386; 
evup [2] “1.0; 


determine if need the default outline or not 
if (otflg == 0) defout (); 
loc = 0; 


Bebeolerı255,255,255); 
elear(): 


/* create sky polygon */ 
bgnpolygon (); 


£or (3=0; j<=ntot;3=3+1) 
( 
cont [x]=(alpts([3] [x]); 
cont [y]=(alpts()3] [y]); 
cont [{z]= -1023.0; 
eral =e(icont{y])* =—90.00058)+0. 880) ; 
es Eleva); 
ot (cont); 


/* create first and last vertices 


Fsszlalets[intot]) [x]): 
if(ts <=525) cont[x]=0.0; 
Pewee 2525) cont [x]=1050.0: 


cont [y)=(alpts[ntot] [y]); 
eeneizj. = -=1023.0; 
elevada: 

y3r (eont); 
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a 


= 


if. (cont |x] == 070) 

( 

e3tilevup), 

var (vupl),; 

eat lezup): 

v3f (vupr); e 

if (alpts[Oj(x) != 1050.0) 
{ 
cont [x] 1050.0; 


contly]*=MoaTpes U a7 
cont [z] = -1023.0; 
c3f (cvd); 


v3f (cont); 
) 


else 
( 
CI (Cvupi. 
v3E(VubE)” 
e3£ (Ev up) 
v3 (Vapi: 
if (alpts[0] [x] != 0.0) 
( 
cont (x] =0.0; 
Conty] =H alpe sko NVD 
cont [zlh=3=1023.0; 
c3£(cvd)” 
VO E (conte) 
) 
) 
endpolygon () ; 


/* initialize the foreground for the mountain range */ 


£frmur [0]=1050.0; 
frmar 1] 150720; 
frmur(2J= -1023.0; 


£frmul [0]=0.0; 
£rmul (1]=40.0; 
frmul(2]= 510230: 


£rmdr(0]=1050.0; 
£{rmdr [1 J=0. 0; 
frmdr [2] =-1023 0; 


frmd1[0]=0.0; 
£frmall 1] 0207 
£rmdl [2 ]e& - 102370: 


£rml(0]=200.0; 


£rml (1]=43.0; 
£frml (2) -2 102320, 
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frm2(0)=400.0; 
Erm2[1)=45.0; 
£rm2[2]= -1023.0; 


pen 03300007 
Eim3 [1)=65.0; 
fem (ic p= -1023,0; . 


£rmá (0]=645.0; 
£rm4(1]=82.0; 
rma {2)= -1023.0; 


£rm5 (0]=790.0; 
ECM =50S. 0; 
Comoras 1023.05 


£rm6[0]=930.0; 
frmbl 2)=130.0; 
frm6(2)= -1023.0; 


rar the change in the z Coordinates for each type of texture */ 


swireh (eurtxtr) 


( 


case 21: /* for rugged granite */ 
RGpEeotor (130,143,128); 

Ze =O 0; 

ZO TA LO. 0; 

zpp = 40.0; 

znn = 210.0; 

break; 


case 22: /* for rocky mountain. */ 
RGBColor (75,42,16); 

zp SO. 0; 

zn 20.0; 

zpp = 60.0; 

zun.= 220.0; 

break; 


case 23: /* £for alpine mountain */ 
RGBcolor (159,88,69); 

zp = 40.0; 

zu =" 25,0; 

zpp = 65.0; 

zune=%230.0; 

break; 


case 24: /* for appalachian mountain */ 
ESBe010r(61,95,36); 

zp 50.0; 

Zhe = 25.0; 

ZEP?- 99.0; 

znn = 240.0; 

break; 


ae, 


case 25: /* fot rolling hal dae 
RGBeEolor (0, 95723017 


Zp = 200 
zn = 930.0; 
zpp = 40.0; 
znn = 190.0; 
break; 


) 


/* create foreground */ 


bgnpolygon(); 
v3£f (Erma ie 
v3£f (frmul), 
v3f (frm). 
v3£(£:m2); 
VIE (ÉS) 
v3r(frm4): 
vatifrmo), 
v3£(frm6); 
v3t(frmur). 


v3flfrmar): 
endpolygon(); 


/* initialize variables for fractal mountain range */ 
yfac = (800.0/squfac)/2.0; 

Il=0; 

IC=0; 

Trage 0, 

ct = 0; 


/* create the first base triangle “7 
P1{x] = 0.0; 
Pl{y] = 0.0; 


Pl{z}] = -1020.0; 
P2[x] = 175.0; 
P2[y] = yfac; 
P2(z} = —s6020; 


Bl) = 3508405 
P3[y] 0:0: 
Ez) ="S3400. 0; 


k=0; 
fg=1; 
ilev = 1; 


/* start creation of fractal meuntainsrange 2, 


Mountain generate (Pl, P27 rs alee 
mntrange (first); 


ct = 1; 
vn=1.0; 
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“determine the maximum number of base triangles to create 
maxtri= (squfac * 12)-1; 


=Or (j=1; j<=maxtri; j=j+1) 


{ 


71-0; : 
IC=0; 
Eg=fgtl; 
/* create triangles 1-5, 7-11 in pattern 2% 
iE (fg 6 6&6 fg != 12) 
( 
PAV eo =.) 2 foc) 7 
EL[IM=>P 2 [y]; 
Eliza 2[z2); 
Bal) P3[x)]:; 
P2[y] = P3[y]+((2*yfac) *k); 


P2[z] = P3[2]+(100.0*k): 


PII21+(350,.0”%yn); 
ey); 


P3[x) 
P3 [y] 


/* change the depth of the z value for each texture */ 

Ao ly > P2[y]) 
( 
P3[z])] = P1(2)+2p; 
) 

else 
( 
P3iz] 
) 


A AS ra 


IEC oes Eg < 12) 
( 
AMES y] < PZ[y]) 

( 

Eo cere ee | z)—zp; 

} 


else 


{ 
Po(z) -= P1(z)-zn; 


} 


if {k >=, 1) k=O: 


Pemereate triangle 6 in the pattern for mountiin range */ 
if (fg == 6) 
( 
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PYl[lxz]- = Pais) 
Pliy) = Esser 
Pile q sees (zl 


P3Ix] = ’Pr2 [x]; 
P3[y] = P2[y]; 
P3[z] = P2[z]; 


p2[x] = Pl[x]; 
P2(y] = Plly]+(yfac*2.0); 
P2[zj== P2[z])+tzpp; 


} 


/* create triangle 12 in the pattern for mountain range */ 
if (fg == 12) 
{ 


Pi(x] = P3(x]; 
P1[y] = P3ly]-(yfac*2.0); 
Pi( zip = rs(2z)=znn, 


TP(x] = P2[(x]; 


TP (y] = P2[(y]:; 
TP[z] = P2[z]; 


P2[x] = P3[x]; 
P2[y] P3[y]; 
P2Zfz)] = Patz. 


P3(x] = TP [x]; 
USA TP [y]; 
P3[z2] =TP[z]; 


znn = 100.0; 
n= 10 
k= k+1; 
fg=0; 
} 
ilev = 1; 
/* send the values to create the fractal */ 


mountain generate (Pl,P2,P3,ilev); 
mntrange (first); 


} 


return; 
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PROGRAM = DRAWMNT.C 


K 


this function draws the mountain range outline 


The left mouse draws a line and adds to accumulative object. 
The middle mouse moves to a new position. : 
The program is exited by popup menu on the right mouse. 


include "g1.h" 
#include "device.h" 


/* 


define the world coordinate boundaries */ 


#define WXMIN 0.0 
#define WXMAX 1050.0 
#define WYMIN 0.0 
#define WYMAX 850.5 
#define WZMIN 1023.0 
#define WZMAX 0.0 


#define MOVETO 11 /* opcodes for our complete picture */ 
#define DRAWTO 21 
#define DELDRAW 31 


#define MAXPTS 1000 /* max number of points */ 

extern float pts{MAXPTS] [3]; /* the xyz coords */ 

extern long totops; /* totál number of points in pts */ 
extern float alpts[MAXPTS] [3]; /* editted xyz coordinates to use for 


outline of mountain range el 
Stern int ntot; /* total number of points in alpts E 
extern int flg; /* flag to determine if new outline or not */ 
entern int otflg; /* determines if default outline is used 27 
extern long op[MAXPTS]; /* the operation to perform ay 
drawmnt () 
{ 

nes single line's location... */ 


float single[2]) [3]; 


/* temp wx,wy coords of line. */ 
float wx,wy,Wwz; 


/* value returns from the event queue */ 
short value,valuex,valuey; 


/* temp variable */ 
float junk, 


long i; T noop temp >/ 


—_ 
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+ 


/* device name returned from the event queue */ 
short temp; 


ant bo ¿“used to exit outline seco, 

float veril. /* array used to create outline 2; 

float ter, epyser, /* variables used to check for duplicate points */ 

intra ue /*loop variables * / 

f; put up some help text as the title * / 

wintitle("Mouse values - Draw / New Location / Erase / Exit-KEYBD e 
line must be 

/* continuous and drawn left to right */ 


/* make button hits on the mouse go to the event queue */ 
qdevice (RIGHTMOUSE) ; 

qdevice (MIDDLEMOUSE) ; 

qdevice (LEFTMOUSE) ; 

qdevice (KEYED) ; 


/* queue mouse dial movements */ 
qdevice (MOUSEX) ; 
qdevice (MOUSEY) ; 


/* only report movements of 5 or more to the queue */ 
noise (MOUSEX, 5); 
noise (MOUSEY,5); 


/* if middle button is hit, report mousex + y */ 
tie (MIDDLEMOUSE, MOUSEX, MOUSEY) ; 


/* if Left button is hit? reper] mousex + yaa. 
tie (LEFTMOUSE, MOUSEX, MOUSEY) ; 


/* queue the REDRAW event (i.e. reshape has been called from the 
window manager) 27 


qdevice (REDRAW) ; 


/* attach the cursor to the mouse */ 
attachcursor (MOUSEX, MOUSEY) ; 


/* turn the Cursor on +7 
eurson(); 
zbuffer (TRUE) ; 


/* determine if new or continued outline mf 
if (flg ==1) 
{ 


wx = pts[totops] [0]; 
wy = pts[totops] [1]; 
wz = -1023.0; 
} 

else 


/* initial location for the cursor */ 
wx = WXMIN + (WXMAX - WXMIN) /2.0; 

wy = WYMIN + (WYMAX - WYMIN)/2.0; 

wz = -1023.0; 


122 


mmode (MSINGLE) ; 
Beonrıg():; 


/* set the mouse to be at the world coordinate spot wx,wy */ 
newmouselocation (wx, wy) ; 


Jeeset up the initial line */ 
single[0] [0] WX; 
single (0) [lL] wy? 
single[(0] (2] = wz; 


single[1] [0] = wx; 
single(1](1] wy; 
single{1][(2] = wz; 


/* put a single moveto into the overall picture */ 
torops.= totops + ]; 
op[totops] = MOVETO; 
pts[totops] [0] = wx; 
pts[totops][1] = wy; 
pes [totops][2]1. = wz; 


EX = 2; 
/* wait for the window manager to exit us...*/ 
while(EX != 10) 


{ 


/* while we have values on the Q, process them */ 
while(qtest() != 0) 
( 


switch (qread (&value) ) 
{ 
case KEYBD: 
if (value == 101). /*e*/ 
( 
/*totop3s = totops =-1;*) 


/* make the single line back one point */ 
single[0] [0] = pts{totops][({0]; 
single(0](1] = pts[totops] (1]; 
single[0] [2] = pts[totops] [2]; 
single[1] [0] = pts[totops] [0]; 
single [IT un, = pesiltoeeps)] (1); 
single[1] [2] = pts[totops] [2]; 


/* set new point to the previous point */ 
op[totops+1]é= op[totops]; 

pts[totops+1] [(O]=pts[totops] [0]; 
pts[totops+1] [1]=pts[totops] [1]; 
pts[totops+1] [2]=pts[totops] [2]; 


EX = 10; 


} 
break; 
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case REDRAW: 


/* should we redraw the picture? */ 
reshapeviewport (); 


/* put the mouse in its new location */ 
newmouselocation (wx,wy) ; 
break; 


case MIDDLEMOUSE: /* if middle MOUSE, start a new line */ 


/* next 2 records on Q are mousex and mousey */ 
temp=qread (&valuex) ; 
temp=qread (&valuey); 


/* compute the world coordinate */ 
computeworldcoord(valuex, valuey, &wx, &wy) ; 


/* make the single line be here */ 
single[0] [0] = wx; 
single(0] [1] = wy; 


single[0] [2] = wz; 
single[1]([0] = wx; 
single{1]{1] = wy; 
single{1] [2] = wz; 


/* put a moveto into the big picture */ 
totops = totops + 1; 
op[totops] = MOVETO; 
ptsítotops] [0] = wx; 
pts[totops] [1] = wy; 
pts[(totops] [2] WZ; 


break; 


case LEFTMOUSE: /* if left MOUSE, draw the line */ 


/* read the queue for the mousex + y */ 
temp=qread (&valuex) ; 
temp=qread (&valuey); 


/* compute the world coordinate of the mouse */ 
computeworldcoord (valuex, valuey, &wx, &wy); 


/* make the single line be here */ 


single(0](0] = wx; single[0] [1] = wy; 
single[0] [2] = wz; 
single(1](0] = wx; 
single[l] [1] = wy; 
single[(1](2] = wz; 


/* put a drawto into the big picture */ 
LOESPSZ=ZZESEOpSE 


op[totops] = DRAWTO; 
pts[totops] [0] = wx; 
pts[totops] [1] = wy; 
pts[totops] [2] = wz; 


break; 
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case RIGHTMOUSE: /* If right mouse delete last point */ 


case MOUSEX: 


case MOUSEY: 


default: 


EGEODS = Cotops =i. 


/* make the single line back one point */ 
single[0] (0) = pts[totops] [0]; 


single(0](1] = pts[totops] [1]; 
single[0][2] = pts[totops] [2]; 
Ssingle[1]([0] = pts[totops] [0]; 


single[1] [1] 
Single([1) [2] 


pts(totops] (1); 
pets (cottons) (2); 


op[totops+1] = op[totops]; 

pts[totops+1] [0]=pts[totops] [0]; 
pts[totops+t1] [1]=pts[totops] [1]; 
pts[totops+t1] [2]=pts[totops] [2]; 


break; 
/* we have movement on the mouse in the X direction */ 


/* get the world coordinates */ 
computeworldcoord (value, valuey, &wx, &junk) ; 


/* must change the x coord of the single line */ 
single[({1])[0] = wx; 
break; 


/* we have movement in the Y direction on the mouse */ 


/* get the world coordinates */ 
computeworldcoord (valuex, value, & junk, &wy) ; 


/* must change the y coord of the single line */ 


single[l] [1] = wy; 
break; 


break; 


} /* end switch statement */ 


) /* end while there is stuff on the Q (gtest() != 0) */ 


/* rebuild the complete picture */ 
ReBeolor(0,175,255); 


clear(); 


/* draw the lines 
linewidth (2); 


thicker */ 


A draw the big picture in black */ 


RGBcolor (0,0,0); 
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/* draw the big = pleeugem a 
for(i=0; 4 <= €cCtops; mesi 
{ 
bgnline(); 
switch (op[i]) 
{ 
case MOVETO: 
endline (); 
bgnline (); 
ver[0]=ptsli] [0]; 


ver[1]=pts[i] [1]; 
ver[2]=ptsli] [2]; 
v3f (ver); 
break; 

case DRAWTO: 
ver (Opes keine, 
ver a) ser 
ver[2]=pts[i] [2]; 
v3f(ver); 
break; 

default: 
break; 

} 
} 


endline(); 


/* draw the single line */ 
RGBcolor (295 0-0) 
bgnline() ; 
for (i=0;1<=1;1=1i+1) 
{ 
ver[O]=single[i] [0]: 
ver[l1]=single[i] [1]; 
ver[2]=single[i] [2]; 
v3f (ver); 


) 


endline (); 


} 


/* create a new array without duplicate point and starting 
with the third point. =, 


j=0; 
a O On 
alpesi Jini] 
alpts[j] [2] 


pts [12M ely 
pest +2 77 
-1023.0; 


for (k=3;k<=totops;k++) 

( 
tpx = abs(pts[xk] (0]-alpts([3] (0]); 
tpy = abs (pts[k] [l1]J-alpts[7][1]); 
tp = tpx; 
if (tpy > tp) toic EP 
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wie (tp >=.5) 
{ 


j=j+1; 
alpts[j] [0] = pts{k] [0]; 
alpes 7) Its eine 
alpts[j] [2] = -1023.0; 
) 
| ; 
meOt = 77; 


ICRC Ot <= 3) otflqg = 0; 


Wines Le ("Fractal Mountains") ; 
zburrez FALSE): 


cri er, 


} 
/* put the mouse where the world coordinate last was */ 


newmouselocation (wx, wy) 


float wx,wy; /* last position in world coords for the mouse */ 
{ 
Mona o rrl91nx or1giny; IO rTaluni Ot the Graphics port */ 
long sizex,sizey; /* size of the graphics port */ 
MORE taltx, inity; /* where the mouse should be in the graphics port */ 


/* re-get the location and size of the graphics port */ 
Betorıgin(£koriginx,&originy); 
getsize (&sizex, &Sizey); 
/* we must reset the mouse cursor to be at the location 

of the last x + y coordinate set */ 
initx = ((sizex - 1)/(WXMAX - WXMIN)) * (wx - WXMIN) + originx; 
inity = ((sizey - 1)/(WYMAX - WYMIN)) * (wy - WYMIN) + originy; 


setvaluator (MOUSEX, initx, 0, XMAXSCREEN) ; 
setvaluator (MOUSEY, inity, 0, YMAXSCREEN) ; 


} 


/* compute the world coordinate corresponding to this mouse position */ 


computeworldcoord (mousex, mousey, wx, wy) 
short mousex,mousey; /* the location of the mouse */ 


float *wx, *wy; /* the returned world coord of the mouse */ 
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long Originx, origin. /* origin of the graphics port */ 


long sizex, sizey; /* size of the graphics porcii 


/* re-get the locätion and size of the graphics port */ 
getorigin(&originx, &originy); 
getsize (&sizex,&sizey); 


/* compute the world coordinate */ 
*wx = ((WXMAX - WXMIN)/(sizex - 1)) * (mousex - originx) + WXMIN; 
kwy ((WYMAX - WYMIN)/(sizey - 1)) * (mousey - originy) + WYMIN; 
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PROGRAM GENDYN.C 


Asis file gendyui-c that generates the fractals */ 
include "libem h" 


include "stdio.h" 
include <math.h> - 


*define x 0 
#define y 1 
#define z 2 


/* Global Structures */ 


/* The variables are the following: 
RAND - array for the random variables 
texture - value for each of the different textures 
11,IC - counters for number of base triangles 
Tri - counts # of triangles created 
tail, first -pointers to address of dynamic array 
SEDE EE zeurrent Structural factor 
nlev - total number of recursion levels 
ct - flag for determining if one or more base triangles created 


=. 
extern double RAND[500]; 
extern float texture; 
exrernınt 11,1C,Tri; 
extern LINK. tail, first; 
extern int sStrtfict; 
extern int nlev,ct; 


/* BEGIN RECURSIVE PROCESS */ 
mountain generate (P1,P2,P3,1lev) 


/* Parameter variables */ 
/* P1,P2,P3 -values of endpoint of the triangles 
ilev - current recursion level */ 


float PLL3),E2(3),2373])]: 
int ley: 


/* Local variables 
I - loop variable 
Emidl,Pmid2, Fmid3 =- midpoints of kxk, yeand z coordinates 
TEMP1,TEMP2, TEMP3 - distance squared between endpoints 
ef x and z coordinates 
TEMP - smallest distance squared value between endpoints 
DIST -~ Square root of TEMP 
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TWO - is equal to 2.0 

DPAR - array for putting endpoints into linked list 

yn - flag to determine if triangle small enough 
PTX,PTZ - variables used to determine locatıon ın RAND 
rhead - value returned from ray_to_list 

PT - value to mod to get location in RAND f 
loc - location in RAND / 


int I: 

float Pmidl(3] Pmid2[3] ,Pmid3[3]; 

float TEMP DESTINO: 

float TEMP1,TEMP2,TEMPO; 

float DPAR[{10]; 

float yn; 

float PTX, PTZ; 

/* *ray to list is a function which returns a pointer to a float 
variable */ 

LINK ray tollist (k 

LINK rhead; 


float PT; 
Ine LOC; 
TWO = 2.0; 
/* find the Distance of x and z coordinates of the two endpoints 
of each edge */ 
TEMPO = (P2[x]-P1[x]))*(P2[x]-Pl[x]) + 


(PZ (ZI SPi 2) ete 2 02) ee ee, 
TEMP = TEMPO; 
TEMP1 = (P3[x]-P2[x]) *(P3[x]-P2[x]))+ 
(ESTANSE LS ZII zZeRZ iz, 
if (TEMP1 > TEMP ) TEMP = TEMPI; 
TEMP2 = (Pl[x]-P3[x))*(Pl[x]-P3[x]) + 
(Pl[z]-B3[z)) ITA 
if (TEMP2 > TEMP ) TEMP = TEMP2; 


DIST = sqrt (TEMP); 


if (strtfct == 0) stxrtfict = 16; 
Il=Il+1; 
if (ct == 0) 
( 
af SI DIST <#stertfet) 


= 1.0; 


for (I=0; I<=2; I++) 


{ 

DPAR[I]=P1(1); 
DPAR[1+3]=P2 (1); 
DPAR[I+6]=P3[I]; 
} 

DPAR[9]=yn; 
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if (Il==1) rhead = ray to List (DEAR); 
if (11>Il)ray next (DPAR); 
return: 
) 
) 


/* for the mountain range if it is not the first base triangle stcp 


iterations at the same level as the first base triangle 2 
if (ct == l) 
( 
if (ilev == nlev) 
( 
yn = 1.0; 


for (I=0; I<=2; I++) 
( 
DPAR[1]=P1(1]; 
DPAR[1+3]=P2(1]; 
DPAR[I+6]=P3 [1]; 
) 
DPAR[9]=yn; 


if (11==1) rhead = ray to list (DPAR); 
Yra TI> l ray next (DPAR); 
return; 
) 
} 


/* find the midpoints of x, y and z coordinates */ 


for (120; IT-2; Iti) /* 0 thru 2 => x,y,z */ 
( 
Pmidl[(I] = (P1[(I] + P2{I]) / Two; 
) 


fora (T0; I<=2: Itt)  /* 0 thru 2 => x,y,z */ 
{ 
Pmid2 [I] = (B2[11 + P3[1)) 4 TWO; 
) 


for (I=0; I<=2; I++) /* 0 thru 2 => x,y,z */ 
( 
Pmid3(I] = (P3(1] + P1[I]) / Two; 
) ; 


/* Displace randomly the midpoints of the y coordinate a 


/* for the edge between first and second endpoints */ 
PTX=(P1[(x]+P2[x]); 
PTZ=(-P1[z]+ -P2[z]); 
PI = PT2 + PTX; 
loc" (long) (PT) % 500L; 


r3 


if (TEMPOS S- 0-0) 
( 


Pmidl [y] = ((TEMPO/S50.0)* (texture * RAND [Lloc1)) AER 
} 
else 
{ 
Pmidl[y] = (texture * RAND(loc]) + Pmidl [y]; 
) 
/* for the edge between second and third endpoints “7 


PTX=(P2[(x]+P3(x]); 
PTZ=ei2-E2(zite Esu2)): 
PT = PTZ + PTX: 
loc =v (long) (PT) 3% 500 
if (TEMP1 < 50.0) 
{ 
Pmid2 [y] 
) 
else 
( 
Pmid2[y] = (texture * RAND[loc]) + Pmid2[y]; 
} 


((TEMP1/50.0)* (texture * RAND[loc])) + Pmid2[y]; 


/* for the edge between third and first endpoints yer! 
PTX=(P3[(x]+P1(x]); 
PTZ=( -P3[z]+ -Plfz)); 
PT = PTZ + PTX; 
LOCATION O BE $% 5005: 
inf, (EEMP2 < 5020) 
{ 


Pmid3{y) = ((TEMP2/50.0)*(texture * RAND[loc])) + Pmid3[y]; 
} 
else 
{ 
Pmid3 [y] = (texture * RAND[loc]) + Pmid3[y]; 
) 
yn =7 00; 


for (I=0; I<=2; I++) 

{ 

DPAR(I]=P1[1]; 

DPAR[I+3]=P2[1I]; 

DPAR[I+6]=P3[1]; 

} 

DPAR [9] =yn; 
if (Il==1) rhead = ray to List (DIAN 
if (Il==1) first=rhead: 


ilev = ilev + 1; 
if (fil) cay were PAR 


192 


/* Recurse on the triangles to create smaller triangles */ 


TEL A TE + 1: 


mountain generace (EMidl, BZ, Pmid2,ilev); 


Drets l: 


mountain generate (P3, Pmid3, Pmid2, ilev); 


Tri ="Tri +]: 


mount ainscenerate(emidl,Pmid3,Pl ‚ılevy), 


Tri =F Tri +1: 


if (ct == 0) nlev = ilev; 


mountain generate (Pmidl,Pmid3, Pmid2,ilev); 


/* END generate */ 
return; 


} 


PROGRAM : RAYLT.C 


/* This file raylt.c dynamicly stores the first set of points 


for creating the mountains. 
It transfers the arrays into a list.*/ 


include Tlist.h" 
#include <math.h> 


extern LINK tail; 
LINK ray_to list (DPAR) 


float DPAR[]; 


( 


LINK head=NULL; 

EINK oLldst ° 

Ine ai 

if (DPAR[(O)>=’\0') /*base case*/ 


( 

head = (LINK) malloc(sizeof (ELEMENT) ) ; 
head ->AP = DPAR[O); 

head ->prev = NULL; 


tail = head; 


/* store the rest of the 9 values of the DPAR values 


LOL (al <a 9: 14+) 
( 
tail -> next = (LINK)malloc (sizeof (ELEMENT)); 
oldst = tail; 
tail = tail -> next; 
tail -> AP =DPAR(i]; 
tail =>) prev = oldst; 


} 


return (head); 


} 


133 


a, 


PROGRAM : RAYNT.C 


/*This file raynt is to dynamicly store the points for creating 
the mountains after the first triangle 15 stored with raylt.c*/ 


#ınelude "liste. ho 
extern LINK tail; 
extern int IC; 


ray next (DPAR) 


float DPAR[]; 
antic ie: 
LINK oldst; 
for (i=0;i <= 9;itt) 
{ 


tail -> next= (LINK)malloc (sizeof (ELEMENT) ) ; 
oldst = tail; 
tail = tail -> next; 
tail -> AP=DPAR[i]; 
tail -> prev = oldst; 
) 
return; 


} 
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PROGRAM : INTERPMNT.C 


/* This routine interprets between the points and then takes care 
of backtracking and intersections */ 


#define x 0 
#define y l 


include "list.h" 
#include "stdio.h" 
#include <math.h> 


/* The variables are the following: 
IC,Il - counters for new base triangles 
alpts - array for values of a given outline 
ntot - total number of points in alpts array 
ER 
extermeirnte -1C, 11; 
externt loat alpts[1000] [3]; 
exeern ıntentot; 


anteremnt () 


{ 
Po eEeocal Variables 
1, k - loop variables 
j - variable to increment from 0 to 1049 
ix - last value of each interpolation set 
kl - variable used for interpolation 
intout - array to hold interpolated values 
sl - slope between any two points 
it - variable used in replacing the largest point into 
array alpts 


2, 


PoE I moc A 

ent ınFr9ur [1050] [2]; 
Sloat si; 

iat cit: 


for (i=0;i<=1049;i=i+1) 
( 


= 


0; 


Iineout[ı])[(x] 
ine out 1} (y} 


) 


etar at the far befthand side */ 
O 
ieee atots (5) (x) ); 


/* if the outline started on the left side use this if statement 


APO O) 
( 
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/* create intout array from the x values of pointe!) to pont Mariara 


alpts */ 
intout [ix] [y] = (int) (alpts[3] [y]); 
intout(31(Yy1.= (int MMalpets (51 o oe 
sl = (alpts[jlly] — (float) Anrou [else ee 
/* interpolate between points 1 and 2 A, 
kl=0; 


for (k=3+1;k<=ix-1;k=k+1) 
( 
k1=k1+1; 
intout [klily)] = (int) slskl) fF zineeutt lv 


} 


/*create the rest of the y values of intout array */ 
tor (18073 -entot aw 

{ 

j=ix; 

ix = (int) (alpts[i+1l][x]); 


/* if the x values of the two points are equal take the largest 
value of y of the two points a7, 
if (ix == 3) 
{ 
if (Gott MS to at iS aa 
( 
intout [ix] (y]=(int) (alpts[1+11 M0 
) 


/* if the values of the second point in x is smaller than the first 
then check the y values and get the largest = 
ir. a) 
( 


/* if the y value of the second is larger than what is already 
there recompute the y values between the x coordinates of 
the second point and the point previous 2 


if ((int)(alpts[iFi] [y] E AAA A 
{ 
intout [ix] [y]= (nt) (aleteo i): 
sl = (alpts[i+1](y]-alpts[i-1] [y]) / (alpts(i+1] (x]-alpts[i-1] [x]); 
j = (int) (alpts[{[i- 1] (21); 


/* do interpolation R 
for (k=j3+1;k<ix;k=k+1) 
{ 
kl=k1+1; 
intout (k] [y] = (int) (s1*k1) + intout(3] [y]; 
} 
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/* if the value is smaller reset the second x Point and go tö the 
next point in alpts ay 


else 
{ 
ix = (int) (alpts[i-1] [x]); 
} 
) 


/* if second of the two points is larger 7 
im (ax > jJ) 

{ 

intout [ix][y] = (int) (alpts[i+1] [y]); 

See spt std alpes 11 1yJ)/(£loat) (1x=3); 


/* do interpolation a, 
kl=0; 
for (k=j+1;k<ix;k=k+1) 
( 
kl=k1+1; 
Icono IRIS (inte) (sl*kl) + antout [3] [y]- 
} 
) 
) 


PomGomecremetnc last point in array alpts to the far right 1049 */ 


J= ix; 

ix = 1049; 

entctoutlız) [x] = ix; 

Arona] = (ant) (alpts[ntot] fy] + 25.0); 


sl = (25.0/ (float) (ix-j)); 


/* do interpolation x 
kl=0; 
for (k=j+1;k<ix;k=k+1) 

{ 

kl=kl+1; 

aneomeik)| [vi = (int) (21*kl) + intout(3] ly]; 

) 
) 


/* replace the values in alpts with the value in intout 
if the value in the array alpts is smaller than 
the value in intout 27 
for (1=0;1<=ntot;ı++r) 
{ 


nte) (alpts[i][x]):; 
Prreenteurtse) [y) > (int) (alpts[i]l[y])) 
{ 
alpts[i][y] = (float) (intout [it] [y])}; 


} 


zerTturrn; 
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PROGRAM : MNTRANGE.C 


/* This routine will pick up the stored values and then 
light the polygons for a mountain range */ 


#include "list.h" = 
#include "stdio.h" 
#include <math.h> 


/* The variables are the following: 
IC,I1 - are counters for number of times subroutines gendyn and listlite 
are used 
frst,tail,ist - pointers for addresses in dynamic storage 
AS(10] - array for storing endpoints of the triangles in 
dynamic storage ad 


extern int ICG pai 
extern LINK frst, tail; 
extern float AS[10]; 


mntrange (head) 
LINK head; 


( 
LINK sc, 


ine: /* loop variables 


IC=IC+1; 
lst = head; 
1f(IC==1) frst = head; 
if (IC==1) lst = tail; 
while (lst != NULL) 
{ 
for (1=9;1>=0;i=i-1) 
( 
AS[i} = lst->AP; 
free (1st); 
1st =7]st- > prev: 
) 


litemnt (); 
) 


return; 


) 
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PROGRAM : DEFOUT.C 
Sirs program initializes the outline for creating a mountain */ 


¿include <math.h> Standar a UNIX include file for math library */ 
weaeclude: “last.h"™ 


inelude “stdio.,h" 


#define x O 
#define y 1 
#define z 2 


/* Global Tables */ 
/* The variables are the following: 
curtxtr - current texture value from the menu 
alpts - array for values of a given outline 
Deere total number of points in alpts array 
counter” tor current point 
a 
extern tnt ceurtxtr; 
extern float alpts{1000] [3]; 
extern int ntot; 


detout () 
( 
/* Local Variables */ 
Tne 1; 


ewaten (cCurtxtr) 

{ 

case 21l: rugged granite Ej 
1=0; 

i = i+l; 
alsts[i] [x] = 

Sir est hlv]).= 250.0; 
1 = itl; 

alpts[i]{x] = 100.0; 
elves y)] = 325.0; 
i = i+l; 

alptsiil[lx] = 260.0; 
Ape [ly] = 550.0; 
i = itl; 

alpts[i] [x] = 330.0; 
alpts{i][{y] = 450.0; 
i = i+l; 

algesıı]) [2] = 502.0; 
alpes(2)][y] = 750.0; 
1 = itl; 

Supe (x) = 675.0; 
aipres/ı])[y] = 345.0; 
i = i+l; 

etoresla lx]. = 900.0; 


| 
© 
oO 
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alpts[i] iy] = 625. 0; 
i = i+l; 
alpts{í] [2] = 1050.0; 
alpts[{i] [y] = 400.0; 
ntot = 8; 
break; 
case 22: /* rocky mountain */ 
1=0; 
1 = itl; 
alpts (1 e = 0.0; 
alptsf{il{y] = 300.0; 
i = itl; 
alpts(i] [x] = 155.0; 
alpts{ility] = 5502-0; 
1 = itl; 
alpts[i] [x] = 210.0; 
alpts[{ij{y]) = 500.0; 


1 = itl; 
alpts({i] (x) = 380.0; 


alpts{i][{y] = 700.0; 
i = itl; 
alptsf[i] [x] = 450.0; 
alpts[i]{y] = 580.0; 
i = itil; 
alpts[i] [x] = 560.0; 
alpts(il[y] = 445.0; 
1 = itl; 
alpts(i]([x] = 746.0; 


alpts({ij ly) = 375.0; 
E 


alpts{i] [x] = 850.0; 
alpts My 1 = 50070, 
i = i+1; 

alpts[i] [x] = 945.0; 
alpts{i]{y])] = 670.0; 


1 = itl; 

alpts[i] =] = 1000.0; 
alpts{i]{y] = 600.0; 
1 = itl; 
alpts[i][x] 
alpesi my] 
ntot = 11; 

break; 


1050.07 
57520» 


case 23: /* alpine mountain x7, 
i=0; 

i = itl; 

alpts [aj] (x) = 0.0; 

alpts[i] iyl = 3500. 

i = i+l; 
alpts[i] [x] 
ape siii 
Is cl 


I l 
D + 
No 
no 
Se, 
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aio] ca 150.0: 

alpts(i][y] = 460.0: 

zu Br 

arpe siil [xz] = 25020; 

aloes illyl = 570.0; 

i = i+l; 2 
alptslillz] = 30070; 

alpes{illy] = 530.0; 

1 = itl; 


apes.) (x) = 395.0> 
alpes (2) fy) := 3275-0; 
i = itl; 
alpts[i] [x] = 465.0; 
eE 395.0; 
i = i+l; 
alprok = 6000; 
ips iy. 455.0; 
i = 1+1; 
aloes (ix) = 725.0: 
empesitiiy) = 575.0: 
i = itl; 
es alos] = 775.0; 
aloes) (yy) = 654.0; 
j = itl; 
alpes [1] 1x] = 990.0; 
See 3750; 
o 
apesi] = 1050.0; 
alpts[i] [y] = 405.0; 
HCOE = 12; 

break; 


case 24: /* applachian mountain */ 
1=0; 

a == ue 

alepts{ij{x] = 0.0; 
crocs ii) = 275.07 
1 = itl; 

alot slills] = 75.0; 
alpesi 305. 0; 
i = 1+1; 

alpes] = 125.0; 
alpesiai iy) = 325.0; 
i = itl; 

apes) (a) = 200.0; 
Soo 355.0; 
1 = itl; 

age) = 275.0; 
alpestzlilvi = 360.0; 
i = i+l; 

ales ia) >= 359.0; 
aistalaly} =.325.0; 
i = itl; 

aler Siili] = 425.0; 
alpts[i][y] = 335.0; 


tal 


i = atl; 
alpts[i]lfxz] = 500T0 
alpestı 12] = 5520: 
i = itl; 
alpts[2]0] 27 =737520; 
alpts[i][y] = 405.0; 
i = itl; 
alpts[i] [x] = 65020; 
alptsiilly] = 42970; 
i = itl; 
alpesli Mkh S 729m0; 
alpts[i] [y] = 420.0; 
1 = itl; 
alpts[i] [x] = 800.0; 
alpts[i] [y] 385.0; 
i = itl; 
alpts [i] (so) = 61530 
alpt say) 365.0: 
ji = itl; 
alpts Wr = 2570: 
alpe SoM rI E AOO 
i = i+l; 
alpts [alle] = 1000.0; 
aletas [1 [Y] = 3Z5.C; 
1= i+l: 
alpts[ı]l[x] = 190950.0; 
alpts[1] [v7 27330202 
ntot — 16: 
break; 


case 25: /* smooth rolling hills x 
i=0; 

i = itl; 
alpts[i] [x] = 
alpts [il [u] = 200.0; 
1 = itl; 

Alpes (“= 75.0, 
alpes a [y] =22570 
i = itl; 

alpts[i] [x] = 125.0; 
alpts[i] [y] = 245.0; 
i = i+l; 

alpts[i] [x] = 200.0; 
alpts[i]I|y] = 26520; 
i = itl; 

alptsti)] [Xx] = 275 2G) 
alpts[i][y] = 274.0; 
i = itl; 

alpts[i][x] = 350.0; 
alpts[i].[y] =.255 300 
i = itl; 

alpts[i] [x] = 425.0; 
alpts[i] [y] = 240.0; 


| 
o 
o 
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i = i+l; 
atpes fi) [21 
alpts[i] [y] 
1 = itl; 
alpts[i) [x] 
alpts[i] [y] 
i = i+l; 
alpts(i][x] 
besa 
i= ihl; 
alpts(i][x] 
alpts(i] [y] 
i= itl; 
alpts[i} [x] 
alpts(i] [y] 
1 = itl; 
alpts[i) [x] 
alpts(i] [y] 
i = itl; 
alpts{1] [x] 
aipts[i] [y] 
i = i+l; 
alpts(i](x] 
alpts[i][(y} 
i = itl; 
sipts[i}] [x] 
alpes (i) iy) 
ntot. = 16; 


break; 


} 


I AM 


500.0: 
230.0; 


=o): 
= 2s Oe 


tot 


650.0; 
Zon aC: 


729280: 


= 265.0; 


800.0; 


= 260.0; 


875.0: 
245.0; 


= 925.0; 
= 230.0; 


/* END defout 


) 


1000.0; 
225.0: 


1050.70; 
235.0; 


dió 
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PROGRAM : LIGHTDEFS.H 

/* This file contains the material, light and lighting model defs 
This is lightdefs.h */ 

/* define the appalachian material light */ 

#define MATFOREST 1 


static float dkgreen[)] = { 


EMISSION, 0.0, 020750707 
AMBIENT, 0.015,.002465 72020252 
DIFFUSE, 0.105, 208sas word, 
SPECULAR, 0.0, 0O ORROD O 
SHININESS, 0.0, 

ALPHA, 1.0, 

LMNULL 


}; 
/* define the granite material light */ 


#define MATGRANITE 2 


Static float ltgray[] = { 
EMISSION, 0.0, 0.0 0.0, 
AMBIENT, 0.400, 0247157, .07505, 
DIFFUSE, 0.405, 02305702200 
SPECULAR, 0.464, 0.364, 0.564, 
SHININESS,207207 
ALPHA, 1.0, 


LMNULL 
le 


/* define the alpine material light */ 


#define MATSUNSET 3 


static float mdorange[] = { 
EMISSION, 2070 EFT 
AMBIENT, 07665; 0.205 20. 00s, 
DIFFUSE, 02725,20240: 7802235, 
SPECULAR, 0.070 30-050 
SHININESS, 0.0, 
ALPHA, O 
LMNULL 
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Ex 


define the smooth rolling hills material anar = / 


+define MATGRASS 4 


static float ltgreen[] = ( 
BIIISSION wee, 0.0, 0.0, 
AMBIENT, AO, 0.905, 0.035, 
PIFFUSE, o 0.994, 0.155, 
SEIECULAR § 0207 0.0, 0:0, 
SATNINASS, 0.0, 
ALPHA, 107 
LMNULL 


y; 
/* define the rocky mountain material light */ 


#define MATROCKS 5 


static float itbrown[] = { 
EMISSION, 0.0, 0.0, 0.0, 
AMBIENT, OS 6S eters 35, 0.035; 
DIFFUSE, 0.585, 0.415, 0.325, 
SEECULAR I ©0.077 07.05. 0.0, 
SHININES>, 0.07 
ALPHA, O 
LMNULL 


7 
/* define the light source for morning*/ 
#define LSOURCM 6 


Sstatiesfloat ightsourcm() = { 
AMBIENT, 0.3, 0.2, 
LCOLOR, 1.00.8, ; 
POSITION Te ToO nS, 0.5, 0.0, 
LMNULL 
ye 


/* define the light source for noon*/ 
#define LSOURCN 12 


Staticwriloae tohtsourcn |) = { 
AMBIENT, 0.6, 
LCOLOR, O 
POSITION, 0.0 
LMNULL 


y 
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/* define the light source for evening*/ 


t!define LSOURCELZT> 


static float lghtsource[] = { 
AMBIENT, 0.4, 
al LCOLOR, IeVr 
POSITION, MiS 
LMNULL 
y; 
/* define the light model for morning*/ 
#define LMODM 7 
static float lghtmodn[] = { 
AMBIENT, QZ 
LOCALVIEWER, 0.0, 
ATTENUATION, 1.0, 
LMNULL 
y; 
/* define the light model for noon*/ 
#define LMODN 10 
static float lghtmodn[] = { 
AMBIENT, 09, 
LOCALVIEWER, 0.0, 
ATTENUATION, 1.0, 
LMNULL 
y; 
/* define the light model for evening*/ 
#define LMODE 11 
static float lghtmode[] = { 
AMBIENT, 023, 
LOCALVIEWER, 0.0, 
ATTENUATION, 2120, 
LMNULL 
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ERSGRAM?: LIGHTINIT.C 


Anas, tilde 139 lightinit.c 
This file initializes the light material, source and scene light 
TOL ISFM 
aAmSehe definitions are in lightdefs.h . = 


include “gl oh" 
#include "lightdefs.h" 


/* function to initialize the mountain material x 


initialmat () 
{ 


/* Make defs for the mountain material x] 


ilmdef (DEFMATERIAL, MATFOREST, 21, dkgreen) ; 
lmdef (DEFMATERIAL,MATGRANITE, 21, ltgray) ; 
imdef (DEFMATERIAL, MATSUNSET, 21,mdorange) ; 
lmdef (DEFMATERIAL, MATGRASS, 21, ltgreen) ; 
lmdef (DEFMATERIAL, MATROCKS, 21, 1tbrown) ; 


} 
a ztwmetiıon tO initialize the light source 2) 


initialls() 
{ 
Make def for the light source */ 


imdef (DEFLIGHT, LSOURCM, 14, 1lghtsourcm) ; 
lmdef (DEFLIGHT, LSOURCN, 14, lghtsourcn); 
Iimdet\DELEICHT, LSOURCE, 14, Lahtsource); 


) 


cameron teazsimitialize the,scene light. */ 
o 


( 
/* Make def for the scene light oy 


lmdef (DEFLMODEL, LMODM, 10, lghtmodm) ; 


lmdef (DEFLMODEL, LMODN, 10, lghtmodn) ; 
lmdef (DEFLMODEL, LMODE, 10, lghtmode) ; 
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PROGRAM : LIGHTMNT.C 


/* This routine will pick up the stored values and then 
light the polygons */ 


#define x 0 
#define y 1 
fdefine z 2 

sS 


#define abs (x) (x)>0 ? (x) : -(x) 


Fine lude gag! in: 
#include "device.h" 
#include "list.h" 
#include "stdio.h" 
#include "lightdefs.h" 
#include "math.h" 


/* The variables are the following: 


IC - counter for current base triangle 
frst,tail - pointers to addresses of the array for 
storing the vertices 
AS[10] - array for storing endpoints of the triangles in 
dynamic storage i 
EN. 
extern int IC; 
extern LINK frst,tail; 
extern float AS[10]; 
litemnt () 
{ 
/* Local variables are the following: 
i - loop variable 
PL1,PL2,PL3 - temporary values used for lighting the triangles 
ysno - flag to determine if triangle small enough 


NL1, NLS, NS - variables to determine the normals of each 
triangle surface 
a 
ne 3, 
/*float temp[3];*/ 
float PL1[3],PL2[37, PRE3 i 


float ysno; 
£tloataNnLi [3] ,NLSP>3] 718 


Ge es 
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/* get the values of the endpoints from the dynamic array */ 


rer (i=0:j<=2:j=i+l) 
{ 


PL1{i] = AS[i]; 
PL2{i] = AS[i+3]; E 
PL3[i] = AS[i+6]; 


} 
ysno = AS[9]; 


/* create the vector normal to the triangle surface SL 

NES wee tet yar ol iy] )*(Pus(z)—-PL1{z]))—( (PL3 ly) sEariy)y (ene [z)-Puil[z))); 
NLS (1]=((PL2[z]-PL1[z]) * (PL3[x]-PL1[x]))-((PL3[z]-PL1[z]) * (PL2(x]-PL1[x])); 
NLS [2] =( (PL2[x]-PL1 [x] ) * (PL3 [y] -PL1[y]))-((PL3(x]-PL1[(x])*(PL2{y]-PL1[y])); 


NS 
NS 


(NLS[0] *NLS[0])+(NLS(1) *NLS(1])4+(NLS[(2] *NLS[2]); 
sqrt (NS); 


NL1(0] 
NL1 [1] 
NL1 [2] 


NLS[0]/NS; 
NLS[1]/NS; 
NLS[2])/NS; 


if (ysno >= 1.000000) 
( 


/* light the triangle and send to the screen */ 
bgnpolygon () ; 

n3fiNL]l); 

v3£(PL1); 

v3f(PL2); 

v3f (PL3); 

endpolygon(); 


} 


return, 
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PROGRAM : TUTSUB.C 


= subroutine that does the tutorial aay 
#define EXIT 95 


tinelude "ol.h® 
#include "device.h" 
#include "list.h" 
#include "math.h" 
#include "stdio.h" 


CUCOTI) 


{ 
/* The local variables are the following: 


tutmenu - integer variable for tutorial menu 

pupval - integer value for pop-up menus 

value - value for reading from the queue 

=: 

int tutmenu; 
int pupval; 
short value; 
/* Create the tutorial menu x7 


tutmenu = defpup("Tutorial %t |On Texture %x11"); 
addtopup(tutmenu, "On Lighting Control 3x 2 75 
addtopup(tutmenu, "On Structural Factor %x13"); 
addtopup (tutmenu, "On Create a Mountain %x14 "); 
addtopup (tutmenu, "On Outline %x15 "); 

addtopup (tutmenu, "On Create a Mountain Range %x16 "); 
addtopup (tutmenu, "Clear %x17"); 

addtopup (tutmenu, "exit $%x95"); 


/* clear the screen with a blue background 27 
R6Beolert(Var7 75,235), 
clear (); 

/* set the color to black for the lettering #7 


RGBcolor (07070): 


/* name the title for the window and write out initial introduction 
to the tutorial 
u, 
waneıt leseitrutorıial®e)r- 
cmow214275, 238 
charstr ("Welcome to the Tutorial section for ISFM'"); 
eMOV211(175.650): 
charstr(" This tutorial section will describe the options and defaults"); 
cmovz2i (17956029); 
charstr ("available in ISFM"); 
EMS LL (ISO O 
charstr (" Press Menubutton to see the Contents"); 
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/* stay in loop until EXIT is chosen 2% 
while (pupval != EXIT) 
{ 
switch (qread (&value)) 
{ 


case REDRAW: E 
reshapeviewport () ; 
break; 


case MENUBUTTON: 

if (value == 1) 

{ 
pupval = dopup (tutmenu); 
switch (pupval) 
( 

case 11: /* Tutorial on Texture a 
wintitle("Tutorial on Texture"); 
BGBCOLOr (O 175,255); 
clear (); 
RGBcolor (0,0,0); 
emovzi(275, 800); 
charstr ("Tutorial on Texture"); 
emov21 (495,725); 
charstr(" The texture of the Mountain is the roughness or smoothness"); 
emov21 (1757 1000: 
charstr("of the surface. The possible types of mountains ranging"); 
cmovz2i (1757675): 
charstr ("from the most rugged to the smoothest are: "); 
cemov 2i (22597525); 
charstr(" a) rugged granite with a white and gray color "); 
emev21(225,600)7; 
charstr(" b) rocky mountains with a spiked brown color "); 
cmovz2i (225,575); 
charstr(" c) alpine mountain with a rust and brown color "); 
Emowz22 (225, 550) 
charstr(" d) appalachian with a forest green color"); 
emovZ2i(225, 525) - 
charstr(" e) smooth rolling hills with a grass green color "); 
cmo 216 150,375); 
charstr(" The default value is an alpine mountain"); 
break; 


case 12: 72 Tutorsaleconaricheine®eontrol */ 
ROBe819r.1(0,175, 255); 
clear (); 
RGBeolosr(0,0,0); 
emey21(275,800):; 
cDacsto( Tutorial on Lighting Control"); 
emov21 (175,725) ; 
charstr(” There are three positions the sun can be located"); 
emow21(225,675); 
charstr("a) Morning The sun is located at an angle 30 degrees"); 
emovZ1 (250-7650) ; 
charstr("from the horizontal negative x axis and 30 degrees "); 


Sl 


EMO AO al 

charstr ("out from the xy plane "); 

emov 2225575) 

charstr("b) Noon The sun is located at an angle 60 degrees "); 
Emov2 (2905500 

charstr ("from the horizontal and directly behind the person. "); 
emoav21(225,500%7 

charstr("c) Evening The sun is located at an angle 30 degrees "); 
cmoyv21 (290,475); 

charstr ("from the horizontal positive x axis and 30 degrees "); 
Emovzi1(250, 450)" 

charstr ("out from the xy plane "); 

emoyzaalla0, >79), 

charstr(" The default value is a morning location"); 


break; 

case 13: /* Tutorial on Structural Factor */ 
RGBcoloro, 175,2995 
clear (); 


RGBcolor (0,0,0); 

cmov 2i (27957800), 

charstr ("Tutorial on Structural Factor"); 

cmovz2i(l75, 725): 

charstr(" The structural factor determines how detailed the"); 
cmov2i(176, 700); 

charstr(" surface structure is. The two options are:"); 
cmov2ł1 (225,650); 

charstr("a) Medium The number of recursion levels stop when "); 
cmov2 2509625); 

charstr("the distance of all three edges is less than or "); 
cMOV2 11 2507600) ; 

charstr ("equal to 16 pixels "); 

cmov2i (225,950; 

charstr("a) Fine The number of recursion levels stop when "); 
cMOV2ER 250 752 IE 

charstr ("the distance of all three edges is less than or "); 
emov2a4 2509500): 

charstr ("equal to 9 pixels "); 

cmov2i(150,4251 

charstr("The default structure factor is medium "); 


break; 

case 14: /* Tutorial on create a Mountain =, 
RGBeeo1l8r40,175, 255); 
clear () ; 


RGBcolor(0,0,0); 
cmoy 21 (275,800) 
charstr("Tutorial on Create a Mountain"); 
cmov21 (175-1231 


charstr(" This option is used to create a single mountain."); 
cmov2i (1757700); 


charstr(" The Create a Mountain option allows one to view the"); 
cmoyvy21 (175,679); 


charstr ("different available types of fractal mountains quickly"); 


132 


cemovz1 (290,625); 

charstr(" To execute the create a mountain there is a submenu "); 
emoves (200,600) > 

Slagseh (G2 ving thewmagMiticatton factor, | Fae magnification of "); 
Emovs 1 (2007.57 5) > 

charstr("1.0, 2.0, 3.0 and 4.0 allows one to see a closer view "); 
emEy21%200,550): 3 

GWarstr( Of what the structure really looks like "); 

emovzi (150, 425)" 

eharcst. (Tine defaultzmagnifreation is 1.0 "); 

break; 


case 15: /* Tutorial on Outline 7 
REBeoLlor (175,255), 
elear(); 
RGBColor (0,0,0)? 
emov2i (2757600); 
enarstr( IWiutorialson Outline") ; 
emovci (175,725); 
char str" The outline section is for creating the contour of the.”); 
ENS m2 MO, 100); 
ehairstr ("mountain range. The outline options are to either "); 
emoV23.(175,675): 
Chħharstr ("create the contour of the mountain range or to clear"); 
cemovz221 (175,650); 
charstr ("and reset the contour already created."); 
emovz21 (2007600); 
charstr(" a) To create the outline the following things can be done "); 
emoeu21:(250,5979); 
charstr(" 1) press the left mouse button to draw the outline "); 
emovzı (250,550: 
charstr(" 2) press the middle mouse button for a new location "); 
emov21 (250,525); 
Chħharstz (7 3) press the right mouse button to erase thei points one "); 
cmovy2i (240,500): 
eae Sea at a Came @starting with the last point and working back."); 
emov2r (250,475); 
eHärser(Aä)press the keyboard ’e’ value to exit the outline section"); 
Emoy211X200,450):; 
charstr(" b) Clear and reset- clears the screen display and reset "); 
emoy21 (250,425): 
narco (ano tralizestthe point array to 0.0, so a new outline can "); 
SMmove az 50,400) 
charstr(“ be created. "); 
emov2.(200, 375): 
coe er slo change Ene current Outline go back into create and the"); 
emeyzzi (200,350) ; 
Chħharstr(" cursor will start at the last point created of the previous"); 
evli 2007; 325); 
ebarser (0 outline. "); 
cmovz2i (150,250)? 
Sgarser( default eontour for each type of mountain is available "); 
break; 
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case 16: /* Tutorial on create a Mountain Range */ 
RGBeoalor (0, 1 zes m 

clear (); 

ReBesleor (0,0798 

Emov 21 (275373800) 

charstr ("Tutorial on Create a Mountain Range"); 

cmovzı (L9, 725) 

charstr(" The create a mountain range option is used to create."); 
cmov21 (175 100) 

charstr(" a whole mountain scene."); 

cmov2i(175, 675): 

charstr ("If the outline or the other options are not selected the"); 
cmov224175,620):; 

charstr ("default options will be used. "); 
break; 


case 17: /* for clear */ 
wintitle (“Tutorial”): 
RGBcetler (07 17 S255 )5- 
clear (); 
break; 
} 
} 
break; 
default: 
break; 


} 
} 


/* reset window title and clear the screen */ 
wintitle("Fractal Mountain"); 
RGBcolor (0; 175,2955); 
clear (); 


return, 


) 
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